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1. Introduction 

Understanding the evolution in time of macroscopic quantities such as pressure or temperature 
is a central task in non-equilibrium statistical mechanics. We study this problem rigorously for a 
model of mass-bearing Brownian particles that are prone to coagulate when they are close, where 
the macroscopic quantity in this case is the density of particles of a given mass. Brownian motion 
arises in the particles of a colloid, due to the random agitation of the much smaller molecules 
that form the ambient environment. As such, our model could be considered as one of a colloid, 
where the dominant interaction between particles is that of coagulation. A theoretical discussion 
of coagulation in colloids was undertaken by Smoluchowski in [7]. 

In the model that we consider, a large number of particles, each carrying some integer-valued 
mass, are, at some initial time, scattered in W^, whose dimension d satisfies d > 3 for the purposes 
of this paper. These particles then perform Brownian motions. There is an A-dependent parameter 
e that specifies the range of interaction of any particle in the model: a pair of particles is liable 
to coagulate (to form a new particle that combines the mass of the old two) when the distance 
between the two is of order e. The choice of e as a function of A is dictated by insisting that the 
so-called mean free path is bounded away from zero and infinity in the limit A — > oo of high particle 
number. (The mean free path is the mean time until the first collision of a particle drawn uniformly 
at random at the initial time. A scaling that produces a bounded mean free path is called a kinetic 
limit). Our model incorporates a significant degree of physical realism, absent from earlier work 
on this type of problem, in the sense that we permit the diffusion rate of the particle to depend on 
its mass, including the case where this rate is taken to be decreasing in the mass (it is physically 
reasonable to suppose that the diffusion rate of a Brownian particle is inversely proportional to the 
mass). As we will later describe in precise terms, we also introduce a parameter into the mechanism 
of reaction which allows us to study such reactions over a natural range of their strengths. 

We study the macroscopic evolution of this particle system by measuring the density of particles 
of a given mass m in the vicinity of a macroscopic location x and at some time t. We will prove 
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that, when the initial number of particles is chosen to be high, this density typically evolves as the 
solution of the Smoluchowski system of PDE, 

(1.1) ^(x, t) = d{n)AUx, t) + QUf)ix, t) - Q^{f){x, t) n = 1, 2, . . . 

with initial data /n(-,0) = hn{-), to be specified in more detail shortly. The first term on the 
right-hand-side of (|1.H) corresponds to the diffusion among particles of mass n, with d{n) being 
one-half of the diffusion rate of such particles. The terms in (jl.lj) corresponding to the interaction 
of pairs of particles are given by the gain term 

n 

(1.2) Q?(/)(x,t) = J^/3(m,n-m)/„(x,t)/„„„(x,t), 

m=l 

and the loss term 

oo 

(1.3) Q^if) = 2/„(x, t) J2 Pirn, n)fm{x, t). 

m=l 

Here, the collection of constants /? : — > (0, oo) quantify the macroscopic propensity of mass at a 
pair of values to combine. As well as deriving the system (|1.1|) as the typical macroscopic profile of 
our random model, we prove in this paper the precise relation between the macroscopic constants 
/3 and the microscopic mechanism of reaction. 

We now give a precise definition of the microscopic process. We in fact define a sequence of 
such models, indexed according to the initial number N of particles in them. We define a range of 
interaction, e, according to Ne'^^'^ = Z, where the exact value of the positive constant Z will shortly 
be given. (We will explain heuristically at the beginning of Section |21 why this relation between 
and e determines the regime of bounded mean free path). In defining the model, the main elements 
to describe are the initial random choice of particle locations and masses, the diffusive dynamics, 
and the mechanism for coagulation. 

To describe each of these, we require notation for labelling the particles in this time-dependent 
model. Let a countable set / of symbols be given. A configuration q of particles is an x N- valued 
function on a finite subset Ig of /. For any i E Iq, q{i) may be written as {xi,mi). The particle 
labelled by i has mass and location Xj. In practice, the index set Iq will be a function of time, 
with a change occurring only at collision events, of which there are finitely many in any given 
sample of one of the random models. 

As for the dynamics of the process, the action on F of the infinitesmal generator L is given by 

(1.4) {LF){q)=AoF{q)+AcF{q), 

where F : {U.'^ x N}^ [0, oo) denotes a smooth function, its domain being given the product 
topology. In (|1.4j) . the diffusion operator Aq is given by 



(1.5) 



AoF(g) = Y,d{mi)A^^F, 

i&Iq 
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-F{Sy)-F{q) . 
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• the cohection of constants a : — > [0, oo) are the parameters of strength of interaction 
between pairs of particles of given integer mass, to which we earlier alluded. 

• the function V : M.'^ ^ [0, oo) is assumed to be continuous, of compact support, and with 
J^d V{x)dx = 1. Its role is to include among the models we consider a rule for coagulation 
time that may be rather arbitrary, beyond the insistence that it be Markovian and cause 
reaction of a pair of particles at some time when this pair are to be found within an order 
of the range of interaction e. 

• we denote by Sjjq that configuration formed from q by removing the indices i and j from 
Iq, and adding a new index from / to which Sjjq assigns the value (xi,mi + nij). The 
configuration Sf^q is defined in the same way, except that it assigns the value {xj,mi + mj) 
to the new index. The specifics of the collision event then are that the new particle appears 
in one of the locations of the two particles being removed, with the choice being made 
randomly with weights proportional to the mass of the two colliding particles. 

While this choice of new particle location causes a pleasing cancellation of terms at numerous times 
in our proof, it is inessential: the same of method of proof yields our main theorem for any choice 
of location for a new particle that is microscopically close - that is, within an order of e - from the 
place where the two old particles reacted. 

We explain further our reasons for choosing the form of the collision term in (|1.6|) . and the 
interpretation of the various terms. Suppose that two particles (xj,mj) and {xj,mj) are such that 
Xi—Xj enters the support of Ve{-) := at time to- They are likely to remain at a displacement 

of order e for a time period of order e^. Given that the dimension d is at least three, it is likely 
that after a time that is large compared to e"^, the particles have moved apart to a greater distance, 
and that they will not return to an e-vicinity of each other at any later time. This means that the 
integral 

It = a{mi,mj)Vf[xi — Xj)dt, 

Jt=to 

whose derivative in T is the infinitesimal rate at time T for coagulation between the pair, and 
which increases as T rises, is likely to have gained its eventual value as a function of T when T 
is of the order oi Iq + Ce^. Thus, the normalisation by that appears in H1.6|l ensures that the 
integral It reaches an eventual value that is of unit order. Neglecting the effects of other particles, 
the probability of collision between the pair before time T is equal to 1 — exp | — It}- We are 
thus prescribing a dynamics in which the fraction of instances of pairs of particles of displacement 
of order e that result in a coagulation of that pair is bounded away from and 1, uniformly in 
small values of e. The role of the constants a is to determine whether this fraction is high or 
not. The reason then for introducing the constants a is that they range over the right scale to 



while the collision operator Aq is specified by 



(1.6) AcF{q) 



a{mi,mj) 



TTT-i 



F{Slq) + 
rrii + rrij rm + m 



Here, 
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determine precisely how liable mass-pairs are to react microscopically. The choice of the value 
of these constants transmits to the macroscopic reaction propensities appearing in the interaction 
terms ()1.2|) and H1.3|) of the Smoluchowski PDE. The recipe for determining (3 from a is as follows: 
there exists a solution u = Un,m '■ — > (0, oo) of the equation 



(1.7) 



a{n, m) 
d{n) + d{m) 



V{x) 



l+Ur, 



(x) 



that is unique subject to the decay condition Un,m{x) = 0{\x\'^ as |x| 
/3 : N X N ^ (0, oo) in are then specified by the formula 



oo. The quantities 



P{n,m) = a{n,m) / V{x){l + Un,m{x))dx. 



Our main result is conveniently expressed in terms of the empirical measures on the locations of 
particles of a given mass. For each n G N and t G [0, oo), we write gn{dx,t) for the measure on M*^ 
given by 



gn{dx,t)=e'^ ^ J2 ^x,it){dx)l{mi{t) = n) . 



We also require a mild hypothesis on the diffusion coefficients d : N — > (0, oo). Namely, we suppose 
that there exists a function 7 : — > (0, 00) such that a < 7, with 7 satisfying 



(1.9) n27(ni,n2 + ns) max 1, 



d{n2 + n^)^^'^ [d(n2 + 723)] 2^-1-. , \ , \ 

1 < (n2 + n3)7(ni,n2). 



d{n2) 



d{n2) 



The initial random configuration of N particles is formed by scattering particles of numerous 
masses independently in according to densities that are prescribed for each mass. These densities 
will be chosen as continuous functions : M*^ — > [0,oo),n G N}, and should satisfy some fairly 
weak bounds. To be specific, we insist that 

• k* {-A)-^k G L'^^iM.'^) where k := J27=i and k{x) = k{-x). 

• For every m, Y.'^=id{nY^'^^^^'^l{m,n)hn G L'^c{W^) where hn{x) = f hn{y)\x -y\~'^'^'^/^dy. 

• For every m, Y^^^i d{nY/'^~'^/^'^{m,n)hn G L'^JW'') where hn{x) = J hn{y)\x — yl^'^^^^'^dy. 

We then set Z = Yl'^=i /r^ ^ (0, 00) and choose points in M'^ x N independently according 
to a law whose density at (x,n) is equal to hn{x)/Z. Selecting arbitrarily a set of N symbols 
{ij : j G {1, . . . , N}} from I, we define the initial configuration by insisting that qo{ij) is equal 
to the j-ih. of the randomly chosen members of x N. 

Remark It is not hard to show that our assumptions on the initial data {/in} are satisfied if k is 
bounded, k has a bounded support, d(-) is bounded, and 'y{m,n) < C{m)n for a function C(-). 
Indeed if k is bounded and has a bounded support, then k * (— A)~^A; G and hn, hn G Lf^^ for 
every n. It is worth mentioning that if k belongs to the negative Sobolev Space = W~^''^, 
then k * {—A)^^k G L°°. We now state the main result of this paper. 
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Theorem 1.1. Let d > 3. Let J : Mf^ x [0, oo) ^ M 6e a bounded and continuous test function. 
Then, for each n E N and t € (0, oo), 

(1.10) limsupETv 

where we recall that e is related to N by means of the formula Ne'^~'^ = Z , with the constant Z G 
(0, oo) being given by the expression Z = '^^n&i /r^ kl.Kffl . {fn : M'^x [0, oo) [0, oo), n € N} 

denotes a weak solution to the system of partial differential equations iLl\) . with the collection of 
constants /? : ^ [0, oo) being given by iL8\) . 



J{x,n,t)[gn{dx,t) - fn{x,t)dx) 



Remarks 



Conditions for uniqueness among weak solutions of (jLlj) are provided in Proposition 2.6 of 
|inj . Although these conditions are slightly stronger than those we have prove the macro- 
scopic equations to satisfy in our case, we anticipate that this uniqueness holds. Indeed, by 
assuming (for technical rather than physically significant reasons) that inf^^meN 0({n, m) > 0, 
we may verify the condition appearing in Proposition 2.7 of ^01; hereby establishing that 
any weak solution of (jLlf) conserves the total mass YlmeN''^ Ir^ fn{x,t)dx throughout all 
time t G [0, oo), in the case in question (namely, with the coefficients /? : ^ [0, oo) given 
by p.8|) ). A lack of uniqueness in the solution would arise in the presence of a mechanism 
for coagulation that is sufficiently strong to give rise to an effect of gelation, in which this 
total macroscopic mass drops after some finite time. The third among these comments 
briefly discusses how this behaviour might arise from a variant of our model. 
Included in the space of parameter values that satisfy H1.9() is the case where the diffusion 
rate d is a decreasing function of the mass, and the coagulation propensities a satisfy 
a(n,m) < Cum. In fact for a nonincreasing d{-), the condition (|1.9() is equivalent to saying 
that a{n,m) < C{n)m for a function C{n). Also, if the microscopic coagulation rate a 
is identically constant, then the condition (|1.9|) is equivalent to saying that the function 
(i(n)n^/(^~^'^) is nonincreasing. 

We believe that the derivation that we have undertaken is of particular interest, because 
the technique of proof that we adopt is robust under numerous changes in the details of 
the mechanism of interaction. An alternative model which we may treat is that where the 
mechanism of interaction is hard-core. In this variant, the reaction mechanism is deter- 
ministic, a pair of particles coagulating at the moment when their locations differ by 2e. 
The hard-core model with constant diffusion rates d{-) = 1 was treated in dimension d = 2> 
in [Hj, with the macroscopic reaction propensity (3 being identically given by the capacity 
2cap(i?o(l)) of the unit-ball. Although we have not written the proof in this paper using 
a notation that includes the case of hard-core interaction, the method of proof yields the 
results of |5]. Indeed, we can prove the analogue of Theorem 11.11 for hard-core interaction 
in any dimension d > 3, provided that we may invoke the hypothesis 



(1.11) 722 max < 1, 



-d{n2 + ns)- 


■id-2 
2 


'd{n2 + ns)- 


- d{n2) 


5 


L d{n2) J 



2d-l- 



< 77,2 + na. 
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The relation [5{n,m) = (d(n) + d(m))caf3(-Bo(l)) holds. The hard-core model is anyway 
approximated in terms of its macroscopic evolution by the stochastic one, with high choices 
for the constants a, as we will see in Theorem 16.21 if the support of V is the unit ball 
Bq{1) and n, m G N are given, then the limit limQ,(„^m)^oo /?(^; w-) exists and equals {d{n) + 
(i(m))cap(i?o(l)) ) the limit being taken with d{n) and d{'m) fixed. 

• Our technique of proof also yields a kinetic limit derivation for the model in which particles 
are assumed to have a range of interaction that is mass-dependent. To give an example of 
such variants, suppose that each particle of mass m has a radius r{m), where r(m) = m^, 
with X a fixed value satisfying x € (O, l/{d — 2)) . Assume, for the sake of brevity, that the 
diffusion rate d is uniformly bounded above and below. We stipulate that particles of mass 
m and n are liable to react when their displacement reaches the order of (r(m) -|- r(n))e. 
More precisely, we modify the definition (|1.6|) of the collision operator Ac by replacing the 
appearance of V by (r(n) -|- r(m)) ■ /{r{n) -\- r(m))), (the factor that multiplies V 
being introduced so that, roughly speaking, the altered collision mechanism respects the 
spatial-temporal scaling of Brownian motion). The statement of Theorem 11.11 is altered 
in that the macroscopic reaction propensities are given by (1.8) with Un,m solving (1.7) 
where V is replaced by (r(n) -|- r{m)^ ^y( • /(r{n) + r{m))^. (The condition that x be 
bounded above as stated is required in the proof of Lemma ESI Varying the details of this 
proof is the principal step required in extending our investigation to this altered model.) 
We mention also that investigating this model in the case where the radius r is chosen to 
be a more rapidly increasing function of mass may be highly interesting, because it may 
be used to describe interacting systems of polymers having a fractal structure, in which a 
gelation effect occurs macroscopically, in the sense that the system of equations (|l.lj) lacks a 
well-defined solution after some finite time before which particles of high mass acculumate. 

• In Section |21 we will discuss how the product of densities that appears in the interaction 
terms (|1.2j) and (|l.^-i|) of the macroscopic evolution is a consequence of the asymptotic 
independence for each positive time of particle locations in the large random models. As we 
will discuss, this independence fails to hold in a large microscopic vicinity of the location 
of a given particle: it is less probable that at any given positive time, a particle is present 
at a distance of the order of e from a given particle than at a randomly chosen point 
in space, because the mechanism of coagulation has the effect of clearing space about 
any given particle. The endeavour of deriving the relation (|1.8() between the macroscopic 
and microscopic reaction propensities has been a formidable one, since in so doing, we 
have determined precisely the nature of this negative microscopic correlation that forms 
a correction to the overall asymptotic independence of particle locations at positive times. 
The authors of [Ij accomplished an analogous task for the quantum mechanical problem of a 
Bose-Einstein condensate with a quite singular choice of pair potential, in which the negative 
microscopic correlation corrects the Born approximation. An assumption that removed 
interaction between particles in small regions where several were present was employed. 

We conclude this introduction by mentioning some previous work related to the problem. 

As we just mentioned, the project of deriving the kinetic limit of the particle densities that we 
undertake was treated for the hard-core model, in the case where the diffusion rates are identically 
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equal to a constant, in [Hj. The approach of this paper was that by which Lanford 4^ had derived 
for short times Boltzmann's equation as a similar limit of a system of balls that move at constant 
velocities and experienced elastic collisions. That is, the authors verified that the correlation 
functions of the particle locations satisfied the BBGKY hierarchy, a system of equations that may 
be derived formally for such systems by assuming sufficient independence in the location of particles. 
Similarly to Lanford's approach, the derivation in |^ was made on a short time interval. However, 
an iteration allowed it to be extended for all time. 

A related problem has been studied by Sznitman [H]. Brownian spheres diffuse and annihilate as 
soon as they touch. The partial differential equation by which the density of particles evolves was 
derived for the kinetic limit, in each dimension d>2. 

The stochastic coalescent is a random model in which particles - elements of an abstract set - are 
assigned a non-negative mass lying in [0,oo), and pairs of particles of mass x and y coagulate at 
an exponential rate K{x,y). Norris jH] showed that for many choices of interaction rate K, in the 
limit of large initial particle number, the density of particles evolves in time according to a spatially 
homogeneous analogue of where the densities depend on time alone and the Laplacian term is 
absent. The coagulation rate P in this case is the same as in the random model: P{x, y) = K{x, y). 
The derivation was valid provided that this system of partial differential equations had an unique 
solution, sufficient conditions for which were proved. 

In [2], a derivation of a system similar to is made from a microscopic model with a finite 
number of possible mass types, and a mechanism for interaction that is over a long-range compared 
to the microscopic scale. 

In the next section, we give an overview of the proof of Theorem ll.il describing the organisation 
of the paper as we do so. 



2. A SKETCH OF THE PROOF 

An essential element in the proof of Theorem 11.11 is to verify that the distribution of particles 
of two distinct masses is independent in the high N limit. This is the so-called Stosszahlansatz of 
Boltzmann. Before sketching how we perform this task, we illustrate how assuming a form of the 
Stosszahlansatz tells us that the choice of scaling = e'^~'^Z in the model identifies the regime 
of the bounded mean free path: that is, why the number of collisions experienced per unit time 
by a typical particle is of unit order, when A^ is taken to be very large. In this regard, recall that 
the volume swept out by a ball of radius e about the location of a particle in the first unit of 
time is that of the 'Wiener sausage' of parameter e, and is of the order of e'^~^. The number of 
collisions is of the order of the number of instances that, in this unit of time, other particles lie in 
this swept-out volume. The other particles are distributed independently at later times - at least, 
if we admit some form of the Stosszahlansatz - so that we expect the particle to encounter of the 
order of Ne'^~'^ other particles. The scaling of A^ and e ensures that this remains bounded above 
and away from zero as A^ increases. 
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In the form in which we prove it, the Stosszahlansatz asserts that the total 'propensity to 
coagulate' between particles of mass Mi and M2, weighted spatially by some time-dependent test 
functions J, J : M"^ ^ [0, 00) x N x [0, 00) in the first T units of time, 

(2.1) e''"^ / Ve{xi - Xj)a{mi,mj)J{xi,Mi,t)J{xj,M2,t)li^mi = Mi,mj = M2^dt, 

is close, when N is high enough, to an expression involving a space-time integral of products of the 
candidate densities: 

rT r 

f6, 



(2.2) p{Mi,M2) / / f{Mi,x,t)f{M2,x,t)J{x,Mi,t)J{x,M2,t)dxdt. 

Jo JR'i 

The precise statement of this claim appears in Section [|-{. 51 as Proposition ^ In seeking to verify 
the claim, we introduce the random variables Q{z), or, more fully, Q{z, Mi, M2,t), given by 

(2.3) e'^-^ V,{xi{t) - Xj{t) + z)a{m,{t),mj{t)) 



j{xi{t),Mi,t)j{xj{t),M2,t)l^mi{t) =Mi,mj{t) = M2}, 



where z G W^. In the case of z = 0, (5(0) measures the average propensity per particle at time t for 
particles of mass Mi and M2 to interact (with a spatial weighting due to the test functions). It is an 
average due to the appearance of the factor of e"^"^ in 1)2. 3|) . whose reciprocal is of the order of the 
number of particles in the system initially. The total number of particles is expected to remain of 
this order at any given later time, because a typical particle collides with only finitely many others 
per unit time. This means that, ignoring a constant factor, we may think of (5(0) as the average 
propensity to interact for any given value t of the time parameter. We also consider expressions 
Q, defined by the formula in (|2.3|) . with the interaction kernel being replaced by another, V^j 
having similar properties. If z is fixed at a small and constant value, while the parameter e is taken 
to be much smaller, then any such Q{z) measures the extent of appearances at time t of particles 
of mass M2 lying in a microscopic ball of radius e centred about a small macroscopic displacement 
z from some mass Mi particle. Instances of such pairs of particles are tentative candidates for 
collision after the passing of a short period of macroscopic time. As such, the assertion that, for 
fixed t € (0, cxo), the expression Q{zi) — Q{z2) is small if zi and Z2 are small, captures the near 
independence of the locations of particles of distinct masses that is in essence the Stosszahlansatz 
. Indeed, we will demonstrate that, for a particular variant Q of Q, 

rT 



0, 



(2.4) limE^ / (Q{z)-Q{0))dt 

zlO Jo ^ ' 

the limit being taken in such a way that e/|z| ^ as 2; — > 0. Then, choosing rf to be a sequence 
that approximates the Dirac b function, we find that Jq Q{0)dt is likely to be close to 

(2.5) f f [ Q{zi- Z2W{zi)r]^{z2)dzidz2dt. 

Jo JR-i JR<i 

By substituting the expression in (|2.3j) (evaluated with the choice z = zi — Z2, and with replacing 
V^) for Q{zi — Z2), we find that the expression in 1)2. 5|) takes the form appearing in (|2.2() . that is, a 
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time-averaged product of candidate densities. (The details, which appear in the proof of Proposition 
n are given in 1)3. 74() ). We now turn to sketch how a statement of the form ()2.4() is derived. 

In a first effort to understand the difference in the behaviour of the time integral of the coagulation 
rate Q{0) and its z-displaced counterpart Q{z), we write down a z-dependent functional Sz on 
configuration space for which the action of the diffusion operator Aq on Xz produces, among 
others, the term Q{z). For each pair {n,m) S N^, we define 4>^^:W^—f (0, oo) by the stipulation 
that A(/)^^(x) = a'(^n,rrije~'^V{z/e), where we write 

(2.6) a'{n, m) = a{n, m) / [d{n) + d{m)) . 
We then set 

(2.7) Sz{q) = e2('='-2) ^ <P%j^^M^{xi - xj + z)J{xi, Mi,t)J{xj, M2,t)l{mi = Mi,mj = M2}. 

When the operator Aq acts, we find that the term Q{z) — Q{0) is obtained when the operator Aq 
acts on Sz — So with each of the two derivatives of the Laplacian falling on c/)*^ rather than the test 
functions. Consider the identity 

(2.8) {Sz-So){T) = {Sz-So){o)+ £ i^^^+AoYsz-Somdt 

+ r AciSz - Somdt + Mt, 
Jo 

and note that the process {Mt : T > 0} is a martingale. We know that two terms lie inside this 
identity that remain of unit order as we take a limit in low z (with e tending to much more 
quickly). These are Jq Q{z)dt and — /q Q{0)dt. Were the only terms to remain of unit order to 
be those two, we would deduce a conclusion of the form we seek, namely 1)2.41) where Q is simply 
Q. However, 1)2. 4|) cannot be valid for this choice of Q. Indeed, imagine that a{Mi, M2) is a very 
large constant. Particles of order e apart would typically coagulate in a very small fraction of the 
order of time in which they are likely to remain at a distance of order e. At the moment that 
the pair of particles vanishes, it ceases to contribute to the expression Q(0), meaning that, at any 
particular time t € [0, 00), Q(0) need not be large, even if the constant a{Mi,M2) is. However, for 
fixed z, pairs of particles, one of which lies in an e-ball about the z-displacement of the other, may 
remain in this state for a significant constant multiple of of upcoming time. This means that 
Q{z) is likely to be huge, if a(Mi, M2) is. 

In fact, there is another term that remains of unit order in p. 8)1 in the low z limit. This term 
arises from the action of the collision operator Ac on 5*0. Two or three particles may index the 
summands of such terms. Two such particles are those that may coagulate, appearing in the sum 
associated to the collision operator. There may or may not be another occurring as a third particle 
in the term Sq. It is the case of the term corresponding to only two particles that remains of unit 
order. It takes the form 

(2.9) e2('^-2)a(Mi,M2) ^ V,{xk - xi)cP%j^^M,ixk - xi)Jixk, M^)J{xi, M2). 

k,l£lq 
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This term witnesses the abrupt curtailment of the propensity to coagulate of a pair of particles at 
that moment when they do coagulate. Given the argument just presented against the possibility 
that (5(0) is equal to Q(0) in (|2.4|) . it is not surprising that this extra term should arise in this way. 

There are numerous other terms in ()2.8() . It turns out that each of them vanishes in the limit 
of low z. Accepting for the time being that this is so, we have however not quite arrived at a 
conclusion of the form we seek. What we have found is that J^^ Q(z)dt is close, when z is small, to 
the expression 



(2.10) e2{'^-2)a(Mi,M2) dt V V,{xk 

Jo : 



Xlj 

1 - ^lh,M2 i^i - ^j)] Ji^i^ Mi)J{xj,M2)l{mi = Mi,mj = M2}. 

The '1' that appears in the square bracket corresponds to (5(0), and the other term to the unit 
order term 1)2. 9(1 . It was our aim to show that Jq Q{0)dt was close to an expression of the form 
Jq Q{z)dt, but our first try gave us such a statement with the roles of Q and Q reversed. We are 
led therefore to consider an altered functional of the configuration space: 

(2.11) X,{q) = £^('^-2) u'm^^mM - Xj + z)J{xi,Mi,t)J{xj,M2,t)\{mi = Mi,mj = M2}. 

We have replaced the appearance of (p^. ^ , (x, — Xj + by that of u%^. ^ , (xj — xj + z) to move 
from (|2.7|) to (|2.1H) . Note that the analogue of the expression (|2.10l) in this case is given by 

(2.12) e2('^-2) [(d(m) + ^(Ms)) An^^^^j^^ (x^ - x,) - a(Mi, M2) V;(x, - Xj)u'{x^ - Xj) 

J{xi,Ml)J{xj,M2)^\rni = Mi,mj = M2}. 

We aim to define the function it^j^^ such a way that this term is equal to (5(0). This means 

that u\,i_^ j^j^ should satisfy 

(2.13) (d(m) + d{M2))^u'M,M^) = a{Mi,M2)Ve{x)u\x) + a{Mi,M2)V'{x) 

for each x G M"^. Writing u\.^^ m2^^) ~ ^^~'^^Mi,M2 (^^/e) and recalling the definition ()2.6|) of the 
constants q'(Mi,M2), the equation 1)2. 13|) for any, or all, e > 0, is equivalent to the requirement 
that umi,M2 : IR'^ — > ]R satisfies 

(2.14) AuM„M2(.x) = a'{Mi,M2)V{x) [l + UM„M2ix) 

In Section El we prove that the equation 1)2. 14() has an unique solution in C'^(R'^). This piece of 
theory allows us to define the functional Xz so that it has the desired effect: the term Jq Q{0)dt is 
close to Q{z)dt for small z. The form of Q{z) is given by 
(2.15) 

^2(d-2) Y {d{Mi) + d{M2))Auli^^M^{x, - xj + z)j{x„Mi)J{xj,M2)l{mi{t) = Mi,mj{t) = M2]. 

This expression explains the form of the modification as we pass from the microscopic propensity 
to coagulate, given by a{Mi, M2), to its macroscopic analogue (3{Mi, M2). The recipe (|1.8|) for 
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computing (3 is nothing other than 

/3(Mi,M2) = ((i(Mi) + d(M2)) / AuM„M2{x)dx 

an equation that arises because for small z, the significant behaviour in (|2.15jl is given by that of 
u{x — y) on its diagonal x = y. 

Our task then becomes that of examining the analogue of the identity (|2.8|) in the case where Sz 
is replaced by X^. Estimates are required to prove that the remaining terms, including that of the 
martingale, vanish in the low z- limit. It would suffice in this regard to demonstrate the following 
form of independence in the behaviour of the particles: that there exists a large constant C such 
that, for any n G N and t G (0,oo), if n particles are picked independently at random from the 
initial configuration, then the probability that all of them lie in a given small ball of radius r at 
time t is bounded above by Cr". Although this statement may seem no easier to prove than the 
Stosszahlansatz , there is a special case for which it may be obtained by a coupling argument. This 
case is the one treated by 0, when all of the particles are supposed to have the same diffusion 
rate (although the argument presented here does not require the mechanism of coagulation to be 
of hard-core type). With this assumption, consider an altered model in which there are particles 
of two types, those living and those perished. Each particle is living at the beginning. Collisions 
may occur only between living particles. When a collision occurs, one particle remains alive and 
the other perishes, with each particle having an equal chance of surviving, independently of other 
randomness. Particles of either type diffuse according to independent Brownian motions at one 
given rate. The ensemble of all particles is a system of independently evolving Brownian motions, 
and is at any given time a superset of the living particles, which evolve according to our collision 
dynamics. The estimate mentioned follows straightaway from this fact. 

In general, however, no such approach is available. After labelling the various terms arising from 
(|2.8j) at the beginning of Section IHl we instead establish various weakened statements in Subsection 
13.11 of which that in Lemma 13.31 is an example whose proof is now discussed. It asserts that 

c-T 

(2.16) Eat / (it ^ Kn^^n2,ni{xk,xuXi)\[mk = ni,mi = ns) < Ce^~'^'^, 

k,l,iGlq 

where K is an instance of a smooth function mapping M.^'^ to [0, oo), and where C is a constant that 
depends on ni, ns and K. The approach to proving such statements as (|2.16|) has much in common 
with the global one by which (|2.4j) is being obtained. We define a function A^_^ „^ : M^'' — > [0, oo) 
on which the action of the diffusion operator Aq yields the generic summand in (|2.16|1 . That is, we 
insist that 

(d(ni)Ai + d(n2)A2 + (i(n3)A3)^^^ „2 „3(xi,X2,a;3) = -K^^^ ,^^^^^{xi,X2,X3). 

We consider the action of the whole Markov generator Aq + Ac on the sum A''^^ ,^^ ,^^{xi,Xj,Xk) 
over those triples of indices {i,j,k) for which mj = ni and nik = n^. We work inductively in 
the lexicographical ordering in for the pair (ni,n3). The action of Aq on this sum naturally 
provides the expression that we seek to bound in (|2.16j) . The action of the collision operator Ac 
produces three types of terms, that might be called gains, losses or internal terms. The gains 
and losses arise from collisions that produce or destroy particles of mass ni or that contribute 
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to the expression in (|2.16|) in its first or third arguments. The gain term appears in an upper 
bound for this expression, and so must be bounded from above. Such a bound is provided by the 
sum of the loss terms from particles of lower mass, terms which are themselves bounded by the 
inductive hypothesis. The internal term arises from coagulations that remove two particles that 
contribute terms in the left-hand-side of ()2.16|) by means of the unrestricted second argument of 
and produce a new particle that itself contributes such a term. We argue that the internal term is 
non-positive, because the term corresponding to the new particle is less than the sum of the terms 
corresponding to the two that are disappearing. It is at this point that the hypothesis H1.9() on the 
model's parameters is used. 

Some additional lemmas presenting straightforward bounds on the function u*" are given in Sub- 
section \A.2\ For the analogue of (|2.8j) for the process X^, — Xq^ estimates showing that numerous 
terms are negligible in the low e limit are given in Subsection 13.31 The martingale term is similarly 
shown to be small in Subsection 13.41 In this way, the statement (|2.4|) is derived. The derivation of 
the Stosszahlansatz , as described after (|2.4|) . is given in Subsection 13.51 

We aim to show that the macroscopic mass densities evolve according to the system of partial 
differential equations (jl.lf) as follows. In Section |S1 we define the functional 

where J : M'' x N x [0, oo) M denotes a test function. Consider the identity 

(2.17) y,(T) = y,(0) + j\c{Yg){t)dt + +^o) {yci){t)dt + Mt, 

with the process {Mt : T > 0} being a martingale. After demonstrating that the martingale term 
is negligible by arguments similar to those employed for the martingale arising from the previous 
functional X^, we consider the identity obtained from 1)2. 17(1 in the limit of low e. We expect to 
derive the system (jl.lj) in weak form, with the function J playing the role of the test function. 
As discussed after (|2.4p . the Stosszahlansatz enables us to replace the term involving the collision 
operator in 1)1.1(1 by a time-averaged product of microscopic candidates for the density, that take 
the form of the empirical measure for particles of a given mass mollified on some scale 5. To pass 
to the limit in low 5 after e has vanished and in this way obtain the weak form of (jl.l|) . we make 
use of the fact that these empirical measures concentrate in the limit of low e on measures having 
densities that are uniformly bounded with respect to Lebesgue measure. This fact is proved in 
Section El 

We conclude this overview by elaborating the analogy between this work and that of |2j that was 
mentioned in a remark after Theorem ll.il The authors of p. consider a Bose-Einstein condensate 
with a potential between each pair of particles of the form N'^V (^N {xi — Xj)) , these particles having 
locations Xi and xj in 3-space. It is proved that the densities matrices in the problem satisfy an 
infinite BBGKY hierarchy of equations known as the Gross-Pitaevskii hierarchy. In our context, 
the assumption that the Stosszahlansatz holds in the naive form that Q{z) — Q{0) is small for 
small z would lead to a derivation of the Smoluchowski PDE ()1.1|) with macroscopic reaction 
rates P{n,m) being incorrectly computed as being equal to the microscopic ones a(n, m). In |2|, 



13 



an analogously crude assumption that the two-particle density matrices may be expressed as a 
product of one-particle matrices yields the GP-hierarchy with an incorrect coupling constant equal 
to that which arises by making use of the Born approximation. In each of these contexts, it is 
necessary to understand the nature of the microscopic adjustment to particle independence in the 
vicinity of each particle to determine correctly the revelant macroscopic constants. The form of the 
adjustment is very similar between the two problems. It would be of much interest to determine 
why this is so. 



3. Establishing the Stosszahlansatz 



Recall that, for any given pair (n, m) € N^, u'^ 



existence is ensured by Theorem 16 . 1 1 that lies in C 

Au'^„^[x) = a' [n,m) Ve[x) 



u^n.m ■ ^ ~^ [0' ^) is the unique function whose 
2/'TBrf\ satisfying lima;^oo m(^) — and 



c) +V'{x) 

where the constants a' were defined in ()2.6() . We are using the notations 



Note that ut, 



(3.1) 



l"^ [0, oo) is given by 



VJx) 



Rd \x-y 



(n,m) 



■n,m I g ' ~^ 



dy, 



where cq = co{d) = {d — 2)' 



with LOd denoting the surface area of the unit sphere S'^ ^ . We are 



denoting by Un^m the function obtained by the choice e = 1. We present the conditions on the two 
test functions J, J : M"^ x N x [0, cxd) ^ M that appear in ()2.7() . It suffices to work with functions 
that take non-zero values for only one value in the second argument, such functions measuring the 
presence of particles of a given mass. By a temporary abuse of notation, we write 

J{x,Mi,t) = J{x,t)l{m = Ml} and J{x, Mi,t) = J{x,t)l{m = M2} , 

where on the right-hand-side, J and J denote smooth maps from M'^ x [0, 00) to M of compact 
support. We will suppress the appearance of the f-variable when writing the arguments of J and 
J. 

Numerous terms arise when the operators Aq and Ac act on the expression — Xq (recall that 
the functions of configurations X^, indexed by 2: G M'^, were defined in (2.8) ). We now label these 
terms. Unless stated otherwise, we will adopt a notation whereby all the index labels appearing 
in sums should be taken to be distinct. This includes the case of multiple sums. For example, 
J2kiei Sie/ fi^k,xi,Xi) denotes the sum of the evaluation of the function / over all arguments 
that are triples {xk, xi,Xi) where k, I and i are distinct indices in /. Note also that, unless otherwise 
stated, whenever the symbol appears in a summand, we mean n^j^ j^j^. 

Firstly, we label those terms arising from the action of the diffusion operator. To do so, note 
that, for a time-dependent functional F of the configuration space, this action is given by 



d 
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Thus, we label as follows: 



(J-^ + Ao^ (X, - Xo){q{t)) = Hn + + iJia + Hu + H2 + H3 + H^, 



with 



(3.2) Hii = e^^'* "^^ a{mi,mj) (^Xi — xj + z) — (^Xi — xj) J{xi,mi,t)J{xj,mj,t) 



_^2(d 2) a{mi,mj)V^(^Xi — Xj^u'^i^Xi — Xj)j{xi,mi,t)J{xj,mj,t) 
^2(d-2) a{mi,mj)Vf:{xi — xj + 2;)«'^(a;i — + z)j{xi,mi,t)J{xj,mj,t), 



\^Jt {xi,mi,t) J {xj,mj,t) + J {xi,mi, t)Jt{xj,mj,t) 



along with 



(3.3) = 2e^^'' d{mi)J{xj,mj,t) u%{xi — Xj + z) — ul.{xi — Xj) ■ Jx{xi,mi,t), 

H3 



2g2(d 2) d{mj)J{xi,mi,t) ul.{xi — Xj + z) — u%{xi — Xj) ■ Jx{xj,mj,t), 



and 



(3.4) i?4 = e^^'^-^^ ^u^{xi-Xj + z)-u^{xi-Xj) 



i,3&Iq 



\d(mi)l^xJ{xi,rni, t)J{xj,mj,t) + d(rnj)J{xi, rrii, t)AxJ{xj,mj,t) 



where fx denotes the gradient of /, and • the scalar product. As for those terms arising from 
action of the collision operator, 



AciX, - Xo)iq) = G,(l) + G,{2) - Go(l) - Go(2), 



15 



where ^2(1) is set equal to 



(3.5) 



rrik + mi 



u^^Xk - Xi + z)J{xk,mk + mi,t)J{xi,mi,t) 
+ M^(xi - Xk + z)J{xi,mi,t)J{xk,mk + mi,t) 



+ 



mk+mil 



u^{xi - Xi + z)J{xi,mk + mi,t)J{xi,mi,t) 

+ u''{xi - xi + z)J{xi,mi,t)J{xi,mk + mi,t) 
u^Xk - Xi + z)J{xk,mk,t)J{xi,mi,t) 

+ n^(xi - Xk + z)J{xi,mi,t)J{xk,mk,t) 
u^ixi - Xi + z)J{xi,mi,t)J{xi,mi,t) 

+ u''{xi - xi + z)J{xi,mi,t)J{xi,mi,t) 



and where 



(3.6) ^2(2) = -e^^'^ ^) ^ a{mk,mi)Veixk - xi)u''{xk - XI + z)J{xk,mk,t)J{xi,mi,t). 

k,l€lq 

The expression in ()2.9() is an analogue of Go (2). The terms in ^^(l) arise from the changes in the 
functional when a collision occurs due to the influence of the appearance and disppearance of 
particles on other particles that are not directly involved. Those in Gz{2) are due to the absence 
after collision of the summand in Xz indexed by the colliding particles. 

Note that 



(3.7) 



^12 + Go (2 =0. 



The process { {Xz - Xo) {t) :t> 0} satisfies 



(3.8) 



[Xz - Xo) (r) = {Xz - Xo) (0) + + Ao) {Xz - Xo) {t)dt 

+ r Kc{Xz-Xo){t)dt + M{T), 
Jo 
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with {M{t) : t > 0} being a martingale. By using the labels for the various terms that we just 
introduced, we find from (|3.8|1 by use of (|3.7|1 that 



(3.9) 



Hn{t)dt+ [ Hi3{t) 
Jo 



dt 



< \X, - Xo\{q{T)) + \X, - Xo\{q{0)) 



+ I \Hii\{t)dt+ [ \H2\{t)dt+ [ \H3\{t)dt+ [ \H4\{t)dt 
Jo Jo Jo Jo 

+ r \G,il)-Go{l)\{t)dt+ r \G,{2)\{t)dt + \M{T)\. 
Jo Jo 

Since J is of compact support, we have that Xz{q{T)) = for T sufficiently large. 

We denote by Pat the measure on functions from t G [0, 00) to the configurations determined by 
the process at time t. Its expectation will be denoted Etv- We aim to prove the following estimates: 
for each T > 0, 



(3.10) 



^N\Hii\{t)dt < C\z\^, 



EN\H2\{t)dt < CUH+S 



EN\H3\{t)dt < C\z\'i+\ 



EN\H4\{t)dt < C|z|^, 







^EN\G,{2)\{t)dt<c(^Y \ 
EtvI^z -Xo(0)[ < C\z\. 

Later we apply the limit \z\ ^ in such a way that ^ — > 0. We will also show that, for each 
rG(0,oo), 

(3.11) Ejv[m(T)^1 < Ge 



d-2 



3.1. Lemmas bounding collision propensity. We prove three lemmas: 
Lemma 3.1. For any T S [0, 00), 



e En 



fdtY^ 



a(mi,'mj)Vf[Xi — xj) < Z. 
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Lemma 3.2. Let J : M'^ — > [0, oo) be continuous, and let H : M'* — > [0, oo) be a twice differentiable 
function such that —AH = J. Then we have the following inequalities, 
rT 

Eat / e^^*^"^) ^ J{xi - Xj)mimj{d{mi) + d{mj))dt < g'^^'^^'^^Kn ^ H{xi - Xj)mimj, 
rT 

Eat / £^('^^2) — Xj)a{mi, mj)mimj H {xi — Xj)dt < e^^'^^^^EAr H{xi — Xj)mimj. 

Remark The fact that the dimension d is at least three imphes that H >0. 
Lemma 3.3. Assume that the function 7 : — > (0, 00) satisfies 

(3.12) n27(ni,n2 + n3)max|l, ^^{j^^^ ^ }< ("-2 + "-3)7(77-1, 712). 

There exists a collection of constants C : — > (0,oo), such that, for any smooth function J : 
R^'^ — i- [0,00), and any given ni,n3 G N, 
rT 

(3.13) Eat/ dt ^ 7(mi,mj)Ve(xi - Xj) J(xj, j;fc)IL{mi = ni,mfc = 713} 

< e^^^^'^^Cn.^ng ^ / .„2,n3(^l)^2,a;3)/lni(2;i)/ln2(3;2)^n3(2;3)f^a;i(iX2dX3, 

where, also given e > anrf 722 G N, i/ie function na • ^^'^ ~^ [0' ^) ^■^ defined by 

(3-14) ((i(rii)A^, +d(n2)A^.2 + d(n3)A^3)^^^_„2,n3(2;i,a;2,a;3) 

= -7(ni,n2)V;(xi - X2)J(xi,X3). 

Proof of Lemma 13.11 Let X{t) denote the total number of particles present in the model at 
time t. Then, for any T > 0, 

EArX(r) = EArX(O) + / ^N^QX{t)dt+ f ^N^cX{t)dt. 

Jo Jo 

From AqX = 0, it follows that 

pT I'T 

- I ¥.NAcX{t)dt = Eat / V a{mi,m.j)Ve{xi - xj) < EArX(O) = N. 
The result follows from the fact that N = Ze^~'^. □ 



Proof of Lemma 13.21 Set 

Recall the mechanism of the dynamics at collision: the location of the newly created particle is 
one of the two locations of the colliding particles, with weights proportional to the masses of the 
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incident particles. We see that when Ac acts on Xq, all those terms indexed by pairs of particles 
one of which is not involved in the collision cancel. Thus, 

(3.15) Ac-X = — e^*-'^^^^ V^{xi — Xj)a{mi,mj)mimjH{xi — Xj). 
Note also that 

(3.16) AqX = e^^'^"^) ^ /\H{x^ - Xj)mimj{d{mi) + d{mj)). 

From the non-positivity of AqX, apparent from 1)3. 15(1 . and the non-negativity of X, follows 

-En [ AcX{t)dt-EN [ AoX{t)dt <EnX{0). 
Jo Jo 

Substituting the form for AqX and AqX in ()3.16() yields the result. □ 

Proof of Lemma 13.31 Note that 

(3.17) ^^^„2„3(xi,X2,a;3) 



-3d+2 



'^^^^^^ jR-i J^d jRd\ d{ni) d{n2) ^ d{n^) ) ^("i''^2) 

J{z,y')Ve{z - y)dzdydy', 
where cq was defined right after (3.1). Define 

^ni,n2,m{'l) = X] ^n-,,n2,ni{xuXj,Xk)'^{mi=ni,mj =n2,rnk = n^}, 

i,j,kGlq 

and 

n2GN 

Note that Ao^nLna is given by 

Xi, Xj , Xk ) 

n2€Ni,j,keIq 

l{mi = ni,mj = n2,mk = m} 
= ^ 7(ni,n2)V;(xi - Xj)J(xi,Xfc)l{mi = ni,mj = n2,mfe = 713}. 

n2eNi,j,keIg 

Note also that 

where Ini,n3 is given by 

(3.19) Vejxk - xi)a{mk,mi) 



~^ni,mfe,n3 ('''j' ^fc; -^j) ~ ^ni,m;,n3 (-^M -^'j -^i) ^{^i = ni,mj =77,3}. 
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The term Gni,n3 is set equal to 

Til — 1 



(3.20) 



V^{xk - xi)a{mk,mi)l{mk = Mi,mi = ni - Ml} 

Ml=l k,l&Iq 

rrik 



E E [ 



"2eN ijeiq 



nik + mi 



-A: 



ni,n2,"'3 



H ^ l{mi = 722, = ^3} 



"3-1 



+ ^ V;(xfc - x;)a(mfc,m;)l{mfc = M3,m; = 713 - M3} 

M3 = l fc,/G/g 



E E [ 

n2GN i,jelq 



nik + mi 



-a: 



H ■ A {xi,Xj,xi) l{mi = ni,mj = 722}. 

mk + mi ^' J 

The term Lni,n3 may be written as a sum of -^^ni.na ^^'^ ^ni.na^ where -^^ni.ns is given by 
(3.21) y^ y^ V^{xk - xi)a{mk,mi) 

n2GN fc,/G/g 



[X]^"i."2,n3(^i'^fc'^')l{™'i = ni,mk = 722, 772/ = ^^s} 

ieiq 

^nun2,n3iXk,Xi,Xi)l{mk = 721, 772^ = ^2, 772; = 723} 
+ '^^nun2,n3iXk,Xl,Xi)l{mk = 72l,7?2j = 772, 772^ = n^} 



ieiq 



and where -Z>ni,n3 is given by 



(3.22) 



Vejxk - xi)a{mk,mi) 

n2eN k,l£lq 



[ ^2 ^ni,n2,n3iXi^^j^^k)Hmi = ni, 722j = 772, 772^ = 773} 



+ ^2 ^ni,n2,n3(.Xk,Xi,Xj)l{mk = ni,mi = n2,mj = ns} 

The term L\^ ,^,^ arises due to the loss at collision of particles indexed by k and / in the contribution 
to l^ni,n3 of summands where two of the three arguments of A'' are occupied by Xk and xi. The 
term -Z^^^^^a carries those summands where only one of the pair of particles, Xk and x/, appears in 
the arguments of A''. The case where only one of the pair appears, and does so in the unrestricted 
second argument of A"^, is not included here, but instead in the 'internal' term, Ini,n3- 

By the hypothesis we see from the form (|XT7|) of A^ that Ini,n3 < 0. From the 

lemma will follow provided that we establish the following assertion. 
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Claim There exist a collection of constants C : — > [0, oo), such that, for each pair (ni, ns) G 
and for all smooth J -.W^ —>■ [0, oo), and e, T > 0, 

En 



where 



—Eat / Aoy„i,„3(t)(it < Rni,n3, 

Jo Jo 

-^ni,n3 C'nijnaC ^ ^ ^ / ^ni,n2,n3 

(Xl, a;2, Xs) /l^i (Xi)/l„2 (X2)/In3 ix3)dXidX2dX3. 



712 SN ■ 

We establish the claim by an induction on the parameters 111,11.3 G N. We proceed to the generic 
step, at which the claim may be assumed for parameter values (MijAI^), where Mi < ni, or 
Ml = ni,M3 < n^. This step includes the initial case where rii = = 1. From 

Jo Jo 
and the negativity of the terms — Ejv Ln^^nz{t)dt and E^v Ini,nz{t)dt, follows 

(3.23) -f EM{AoYn,,m){t)dt<ENYn,,nM+EN [ Gn,,n,{t)dt. 

Jo Jo 

Note that 

l-T "1-1 r-T "3-1 ,-T 

Eat / Gni,n3it)dt < Cn^ N Eat / LMi,n3it)dt + C„3 N Eat / Ln^,M3{t)dt, 

Ah=l M3=l 

where here we have used such facts as 

^ni,n2,n3 <max{^|^,^|^} ' Aj_„^^„3 for j G {1, . . . , m - 1}. 

We do not require any hypotheses on the parameters in this last step, since there are only finitely 
many possible values of mass for particles that coagulate to form a new particle of some given mass, 
and we can take care of each of these cases by choosing the value of the constants Cm^nz- From 
H3.23() . the inductive hypotheses, the elementary inequalities 

^i,"2,"3 ^ ™{^|^'^|^} ' ^ni,n2,n3 for j G {!,..., m - 1}, 

^ni,"2,i ^maxj^l^,^^} ' fori G {!,..., ns-l}, 

and the fact that the particles are initially independent, we find that the expression 

t-T 

E^(Aoy„i,„3)(t)dt, 







is bounded above by 

Cni.nge^^^"'^^ ^ / A''^.,^.^^,n3{^1^^2,Xz)hn^{xi)hn2{x2)hn3{x3)dxidx2dX3, 

which is the second part of the claim. To establish the first part, note that 



f EN(,Ln„n,){t)dt<ENYn„nM+EN [ Gn,,n,{t)dt. 

Jo 
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Hence, the term to which the first part of the claim refers satisfies the same upper bound as the 
one that was just bounded. □ 

The following corollary to Lemma 3.2 will be used in Section 5. 

Lemma 3.4. There exists a constant C such that for any T G [0, oo), 

rT 

e'^^^ETv / dt 2~] rnimja{mi,mj)Ve{xi — Xj) <C. 

Jo ■~^T 

Proof On account of Lemma 3.2, it suffices to find a pair of functions (J, H) = (J^ , H'^) and a 
positive constant 6 such that J > 0, — AJ = H, e'^~'^H'^{x) > 6 whenever Ve{x) ^ and that the 
expression 

is bounded in A^. For this, let us pick a nonnegative smooth function W of compact support such 
that W{x) > 1 whenever |x| < 2r/ where r/ is chosen so that the support of V is contained in 
the ball i?^(0). We then set J'^(x) = e~'^W{x/e) and H''{x) = e^~'^a{x / e) where the function a 
solves — Aa = W . More precisely, a{x) = co{d) f W{x — y)\y\'^^'^ dy . By restricting the domain of 
y-integration to the case \y\ < rj, it is not hard to deduce that if |x| < rj, then a{x) > {d — 2)~^rj. 
Also, a{x) < c\x\'^~'^ for a constant c. The proof of this is very similar to the proof of the first part 
of Lemma 3.5 below and omitted. The boundedness of the sequence {Xj^} is now a straightforward 
consequence of our assumption k * (— AA:)~^(0) < oo on the initial density. □ 

3.2. Bounds on functionals of Un^m' We will verify the assertions presented in (j3.1()j) . The 
following lemma provides the bounds on the behaviour of the functions {"U^^ : — > [0,oo) : 
(n,m) G N} and other functions that will be used in this section. Recall that li^^(x) = u^{x) = 
e'^~'^u{x/e) where u satisfies (1.8). We choose the constant Co so that V{x) = whenever |x| > Cq. 
Recall that k = J2n iT-hn- 

Lemma 3.5. There exists a collection of constants C : — > (0, oo) for which the following bounds 
hold. 



forxe R", -n„,„(x) < |^ and \Vun,m{x)\ < 
for X G M'^ satisfying \x\ > max {2\z\ + Coe,2Coe}, 



(3.24) 



(3.25) 



< 



Cn,in\^\ 



and 



X 



< 



letting H = Hn^m ■ [0, oo) be given by 



-AHn,m{x) = <^m(x + 2)l{|xl < p}, 
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we have 



Hn,m{xi - X2)k{xi)k{x2)dxidX2 < Cn,m{p + \z\Y 

for z G W^, let = Hz,n,m) ■ ^ [0, oo) be given by 



-AH,{x) 



l{|x[ < p}. 



Then, 



(3.26) 



j Hz{xi - X2)k{xi)k{x2)dxidx2 < Cn,m{p + \z\) ■ 



• for any positive integers n and m and a smooth function J of compact support, there exists 
a constant Cn,m{J) such that for any given z G U.'^, the function ^^^^„2,n3 • ^^'^ [0>°°) 
defined by 

(3.27) {d{ni)A^^ + d{n2)A^.^ + din3)A^.^)A'^^^^^^^^^ixi,X2,X3) 



satisfies 



(3.28) 



-7("-i, "2X^(^1 -X3 + z)e^'^V[ l(ki -a^s! < p)J{x3) 



j A^^,„2,„3(xi,X2,X3)/l„j(2;i)/l„2(x2)/l„3(x3)(ixi(ix2d2;3 



<Cn,miJ)e''-Hp+\z\) 



Proof Throughout the proof, we write u for the function Un,m and a' for the constant a'{n,m) 
of ()2.6|) . The dependence of the constants on n and m arises from that of a' , and is also omitted. 
To prove the first part of the Lemma, we will show that 

(3.29) < -u(x) < cmin <! Ix'^""' 

Note firstly that 

ti(x) = —cqo' 



< |Vn(x)| < cmin 
V{y) 



{\x\'-',l}. 



\x - y\ 



\d~2 



dy 



By Theorem 6.1, u{x) € [—1,0] for all x G W^. From this we find that 

y{y) 



t(x) < Cqo' 



^dy 



\x - y\ 



\d-2 



If |x| is large, then |x — y| > |x|/2 whenever V{y) ^ 0, so that 



\x -y 



K.d \X 



\d~2 



If 1x1 is small, then 



vjy) , . 
\x - yr J\x-y\<c \x-y 



^dy < pdp = — , 



proving the first inequlity in 1)3. 29() . In the same way we can deduce the second inequality in H3.29|) . 
Indeed 

Vn(x) = coa'(2 - d) [ V{y) ""'^ \u{y) + l] dy. 
^R'* X — y M- 



If |x| is large, then |x — y| > \x\/2, and 

x-y 



\x - yl^ 



whereas, if \x\ is smaU, then 

|Vn(x)| < c 

J \% 

Thus, 



y\<ci 



\x — y\'^ 



I Vn(x)| < Co min |l, '^j. 



and so |Vn(x)| < c\x\^ 

To prove (|3.24l) in the second part of the lemma, note that 



(3.30) 



< C / e-'^yf^ 



1 



1 



\d-2 



\x + z — ^ |x — y| 
\x — y\'^~'^ — \x + z — y\'^^'^ 



dy 



\x + z — y\'^~'^\x — y\'^~'^ 



dy 



Note that 
(3.31) 
and that 
(3.32) 



\x + z — y\ 



id-2 



\d-2 



X — y\ ^ \ \x — y \ + \z 



X — y\ — \x + z — y\ < \x — y\ 



d-2 



\d~2 



\x — y \ — \z\ 



d-2 



The right-hand-sides of 1)3. 31() and 1)3.32(1 take the form ca'^ for some a G \\x - 
y\ + |z|] . Note that if V{y/e) / 0, then 

\x — y\ — \z\ > \x\ — Cq€ — \z\ > 0, 

by our assumption on As a result, we have that \x — y\ + \z\ < 2\x — y\, so that 



X - y 



I d-2 



\x + z — y\ 



I d-2 



< C\z\ \x — y 



|d-3 



Hence, 



'{x + z) -u'{x) < C\z\ / V'{y) 



\x - y\ 



;dy. 



\x + z — y\'^ 2 

From the inequality |x| > max{2|z| -|- Coe,2Coe}, we deduce that |a; -|- z — y| > 
\x-y\> \x\/2. We obtain ((H2H). 

In seeking to prove (|3.25() . note that 

x + z - y x-y _ {x + z - y)\x - y\^ - [x - y)\x + z - y\'^ 



\x + z — y\ 



\x - y\ 



\x + z — y\ \x — y\ 
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Note that, for any a G 



(3.33) 



[a + z)\af -\a + zj^'a 



< a 



so long as \z\ < \a\. Given that 



Vu'{x) = co(2 - d)a' / V'{y) 



x-y 



u 



we may apply ()3.33|) with the choice a = x — y to obtain 



\/u^{x + z) - Vu'(x) 



< C\z\ 



V%y) 



dy, 



\x + z — y\ 



dy. 



From the inequality |x| > max{2|z| + Coe,2Coe}, we deduce that |x + 2; — y| > \x — y\/2 and 
|a; — y| > 1^1/2- We conclude that 

Vu'{x + z)-Vu'{x) < 

as required. 

To prove the third part of the Lemma, note that 

\H{x)\<ci [ \x - yf-'^ly + zf-'^my] < p} 
by the first part of the Lemma. Hence, 



(3.34) 



\H{xi — X2)\k{xi)k{x2)dxidx2 



< 


Cl 


(, 


y + z 






'\y\<p 






C2 


1 


\y + z 






'\y\<p 




< 


C3 


(, 


y + z 






'\y\<p 





2-d 



2~d 



xi — X2 — '^k{xi)k{x2)dxidx2dy 



where we used our assumption A;*(— A) ^k £ Lf^^ for the third inequality. This establishes the 
third part of the Lemma. 

To prove the fourth part of the Lemma, we firstly find a bound for Vu^(x) = e^^'^Vu(^x/ej. 
Using the first part of the Lemma, 



Hz{x)<cl \x -y\^ '^\y + z[ ''\{\y\< p] 



it follows that 
(3.35) 



\Hz{xi — X2)\k{xi)k{x2)dxidx2 



< Cl 


(, 


\y + z 




'\y\<p 




< C2 


i 


\y + z 




'\y\<p 





l-d 



V + '^dy < C3(p+ \z\). 



xi — X2 — 2/p '^k{xi)k{x2)dxidx2dy 
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We have deduced ^2^. 

As for the fifth part of the Lemma, let us write J {a) for IL(|a| < p) and define the quantity / to 



be 



(3.36) co{3d) dxidx2dx3 / 7(ni,n2;, , . ' x ' n 

n^^^nj J^'" ^ ^^^^> ^(^2) ^(^3) 

/irM(a;i)/in2(2;2)/in3(2^3)^^^(^' - v' + z)£~'^V {^-^^ J {z - y)J{y')dz'dydy'. 



We write 



I = co{M) / + y-)j{z' -y')J{y')G{z',y,y')dz'dydy', 

where G{z',y,y') is given by 



3d+2 



7(ni,n2) y (^ ^^^(^^^^^ ^^d(n2^^ ^^^(nsV ) hn^{xi)hn^{x2)hn.j,{xj,)dxidx2dx2,. 
Using the elementary inequality abc < {a? + b'^ + c^)'^/^ -^yg deduce that G{z' ,y,y') is at most 

(3.37) 7(ni,n2)(d(ni)d(n2)d(n3))'^/^-^/^/|xi-z'r'^'/' |x2 - yr'+'/^l^rs - yT'^'^' 

Jll ,712,713 

(xi)/i„2 ix2)hn3 {X3)dxidx2dx3. 
From our asumptions on /i„ we deduce that G £ Lf^^. Hence, 

I<C [ u'{z' -y' + z)£-^v(^^^)j{z' - y')J{y')dz'dydy'. 

Note that, for fixed iz',y') G R^'^, 

/ V(^)dy = e'^. 

Thus, 

(3.38) / < C7e^~2 /■ u^z' - y' + z)J{z' - y')J{y')dz'dy' 

< Ce'^-'^ j j u'{z' -y' + z)J{z' -y')dz'dy' 

< Ce'-' [ dy' [- l^^J(z'-y')dz' 

Jk J \z' -y' + z\'^ ^ 



1 

\a\<p+\z\ lar 



^ '^^'^^'^ / ^ Ce°'~2(p + |z| 



where K QW^ denotes a compact set containing the support J, and where we made use of the first 
part of the Lemma in the third inequality. This is the bound stated in (|3.28)) . □ 

3.3. Estimating the terms. 
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3.3.1. The case of H14. Note that 



(3.39) Eat / Hi4{t)dt < Ce^'^'^-'^^En [ V |u"(j;i - x^- + - ^/"(xi - xj) 

ijmi = Mi,mj = M2}, 



where the constant C depends on the L°° bounds satisifed by J, J and then' time derivatives. Hence 



E 



TV 



Hu{t)dt 







<i^l + i^2, 



where i^i is given by 
rT 



■^0 , 1^ 



U 



{xi — Xj + z) — u^{xi — Xj) l{mi = Ml} l{mj = M2} 



i,jeIq.\Xi-Xj\>p 



and K2 is given by 



g2(d-2)jg^ f (it V n'(xi - Xj + z)- u^{xi - Xj) l{mi = Mi}l{mj = M2}. 

Jo . „v-r 



i,j€lq-\Xi-X.j\<p 

Firstly, we treat Ki. Note that 

C|z|e2(d-2) 



(3.40)Ki < ' Eiv / J] l{mi = Mi}l{mj = M2} 



i,j&Iq:\xi-Xj\>p 



C\z\T 



i,jeIq:\xi-Xj\>p 



.d-2 



2^ "^0 < -d^T' 



9(0) 



where the first inequahty follows from the second part of Lemma 13.51 and the final one from an 
hypothesis on the initial conditions of the system. We now treat the term i^2- By writing, 



(3.41) K2 < Ce2(°'~2)]E^ f 

Jo . .^r 



+ z) + U^[xi - Xj) 

i,j(iIq'-\Xi — Xj\<p 

\{mi = Mi}\{mj = M2}, 

we obtain an expression on the right-hand-side which may be bounded by applying Lemma 13.21 
The relevant estimate is provided by the third part of Lemma We find that K2 < C(^p+ |z|) . 
Hence, 



E 



N 







dt<Ki+K2< + C{p+\z\)\ 



Making the choice p = \z\'i+''- , we find that 

Etv / 
Jo 



dt < C|z[d+i 
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3.3.2. The cases of H2 and H3. The estimate of Ejv Jq \H^{t)\dt is derived in an identical fashion 
to that of Ejv Jq \H2{t)\dt. Picking p £M that satisfies p > max {2\z\ + Cqc, 2Coe}, we write 



/ EN\H2{t)\dt < Ri + R2, 
Jo 



where 
(3.42) 

and 
(3.43) 



rT 

I! = 2e2('^-2)]E^ / d{mi)\j{xj,mj,t)\ 

Jo ■ r-T I 1^ 



R2 



2g2(d~2)jE 



AT 



/ d{mi)\j{ci 

Jo ■ v-r I 1^ 

u|(xi - + z) - u%{xi - Xj) \Jx{xi,mi,t)\dt. 



Firstly, we examine the sum Ri. Recalling that we consider test functions J and J respectively 
supported on particles of mass Mi and M2, 



: rrii 



Mi,mj = M2} 



where the lower bound on p allowed us to apply the second part of Lemma [3.51 Thus, Ri < C\z\/p'^. 
Secondly, we bound the sum R2. Note that 

(3.44) R2 < 2e2('^-2)||J,|| ||J||d(Mi)E^y"'^ 

^ i,jeIq:\xi-Xj\<p 

ul.{xi — Xj + z) — n|(xj — Xj) l{mj = Mi,mj = M2} 



rT 
Jo 



where || ■ || denotes the L°° norm and the constant C depends on the test functions J and J. The 
expression p.44p is written in a form to which Lemma 13.21 mav be applied. Doing so yields 

R2 < Ce2('^-2)E^ ^i^j {^oixi - Xj) + H,{x, - xj)) , 

where the functions {Hz '■ z € M*^} appear in the fourth part of Lemma 13.51 From that result and 
our assumptions on the initial data, we deduce that R2 <C(^p+ Hence, 

/ EN\H2it)\dt<Ri+R2<C^-4 + C{p+\z\), 
Jo P 
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for p < 1. By making the choice p = \z\'i+^ , we find that 



[ EN\H2{t)\dt < C\z\^ . 

Jo 



3.3.3. The case of H4. Note that 

EN\H4it)\dt < (d(Afi)||AJ|| ||J|| +d(M2)j|J|| ||AJ||)e2('^-2) 







(3.45) En di y2 ^^^^ = Mi,mj = M2} u^ixi - Xj + z) - u'(xi - Xj) 

This quantity is bounded as in the case of H14. 

3.3.4. The case ofGz{l) — Go(l). We now estimate the term 

En\G,{1) - Goil)\{t)dt. 



To ease the notation, we do not display the dependece of J and J on the variable t. Note that 

rp g 

(3.46) / Ejv|G,(l)-Go(l)|(t)(it< Va, 

•^0 ^=l 



where 



rT 

(3.47) Di = En dt a{mk,mi)Ve{xk - xi)\j{xk,mk)\ 

•'^ k,l&Iq 

g2(d-2) ^ I J(2; . , rnj) I U" {Xk - Xi + z) - {Xk -Xi) , 

each of the other seven terms on the right-hand-side of 1)3. 46|) differing from Di only in an inessential 
way. Given this, the estimates involved for each of the eight cases are in essence identical, and 
we examine only the case of Di.We write Di = + D^, where we have decomposed the inner 
i-indexed sum according to the respective index sets 

{i & Iq,i ^ k, I, \xi^ — Xi\ > p} and {i £ Iq,i ^ k, I, \xk — Xi\ < p}. 

Here, p is a positive parameter that satisfies the bound p > max {2\z\ + C^e, 2Coe}. By the second 
part of Lemma l3.5| we have that 



r<\y\^d-2 
< ^^^^En 

pa i 



rT 

/ dt ^2 a{mk,mi)Ve{xk - xi), 
JO 



where we have also used the fact that the test functions J and J are each supported on the set of 
particles of respective masses Mi and M2, and the fact that the total number of particles living at 
any given time is bounded above by Ze'^~'^. From the bound on the collision that is provided by 
Lemma l3. II follows 
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To bound the term D^, note that is bounded above by 

(3.48) CEn f dt (^i^k,mi)V,{xk-xi)l{mk = Mi}e^^'^-^^J2 

° k,ieiq ieiq 

l{\xi - Xk\ < /o}IL{mj = M2} u^{xk - Xi + z) - u^{xk - Xi)^\J{xi,mi)\ 
j-T 

< CEn dtY, a{mk,mi)V,{xk - xi)\{mk = Mje^^'^-^) ^ 

° k,l£lq ielq 

l{\xi - Xk\ < p}IL{m.j = M2} {u^ixk -Xi + z)\ + |n^(xfc - Xj)|j | J(xj,mj)|. 

Note that the last expectation is bounded by Lemma 13.31 because by our assumption on a, we can 
find 7 such that a < 7 and 7 satisfies the assumption of Lemma 3.3. The upper bound provided 
by this Lemma in this particular application is computed in the fifth part of Lemma 13.51 We find 
that 1)2 < C{p+\z\f. 

Combining these estimates yields 

Di<D^ + D^< C^TT + C{p + \z\f. 

pa 1 

1 2 

Making the choice p = \z\''-+^ leads to the inequality Di < \z\d-+^ . Since each of the cases of 
{Di : i € {1, . . . , 8}| may be treated by a nearly verbatim proof, we deduce that 

rEjv|G.(l) - G^{l)\{t)dt < C\z\^. 
Jo 

3.3.5. The case ofGz{2). Recall that 

G^(2) = -e^^'^"^) ^ a{mk,mi)Ve{xk - xi)u^{xk - XI + z)J{xk,mk)J{xi,mi). 
k,ieiq 

If A;,Z E Ig satisfy K(xjfc — x;) / 0, then \xk — xi\ < Cqc, and so 

\xk - xi + z\ > \z\ - Coe > \z\/2, 
provided that \z\ > 2Cq£. This implies that 

C C 
u^ixk - xi + z) < ■ -j-^ < I |, 

\xk — xi + zl"-'"^ \z\"-~'^ 

where in the first inequality, we used the first part of Lemma 13.51 Applying this bound, we find 
that 

/ EN\G,{2)\dt < En / y2 a{mk,mi)V,{xk-xi)dt 

■^0 1^1 -^0 k,iei, 



1 



whose right-hand-side is bounded above by Ce'^ l\A'^ according to Lemma mi As a result, 
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3.3.6. The case of¥.i\i\Xz — Xq\. We now turn to E^rlX^ — Xo|. Assume that \z\ > CqE. Using the 
first and second part of Lemma 3.5 we have that 

EivlX, -Xo(0)| < C\z\ j j^^hMMhMM\x-y\^~'^dxdy 
(3.49) +C JJ^^ hM,ix)hMM\x - y\''"^H\^ -y\< 3\z\)dxdy, 

+C jj^^hMAx)hM2iy)\x - y + z\'^''^l{\x - y + z\ < i\z\)dxdy, 

where L is a bounded set that contains the support of J and J. Using our assumptions on the 
initial data hn we obtain the bound C\z\ for the first term and C|zp/^ for the other two terms on 
the right-hand-side. This evidently implies our bound on EatIX^ — Xo|(0) in (|3.1U|1 . 



3.4. The martingale term. This section is devoted to proving the estimate ()3.11|) . Note that 
M,(r) = X,{qiT),T) - X,(g(0),0) - £ + X4q{t),t))dt 



is a martingale which satisfies 
where Ai{q,t) is set equal to 



3 X 
MziTf] = Ve,v / Ai{qit),t)dt, 
^ i=i -^0 



j£/q,mi=A/i 



u^{xi - Xj + z)J{xi,Mi,t)J{xj,M2,t) , 

jeIq,mj=M2 

with A2 being given by the same expression with instances of Mi and M2 being interchanged, 
including those implicitly appearing in the function = u\j_^ j^.j^ . The term A3 is given by 



^4{d 2) a{mi,mj)V,{xi- Xj) 



El 



r rrii 



nii + m 



'{xi - Xk + z)J{xi,mi + mj)J{xk,mk) 



(3.50) 



+ 
+ 
+ 



rrii + 
rrii + n^-j 



rrii + m 



-u^ixk - Xi + z)J{xk,mk)J{xi,mi + mj) 
-u^{xj — Xk + z)J{xj,mi + mj)J{xk,mk) 
-u''{xk - Xj + z)J{xk,mk)J{xj,mi + mj) 



u^ixi - Xk + z)J{xi,mi)J{xk,mk) - u^ixk - Xi + z)J{xk,mk)J{xi,mi) 
u'^{xj — Xk + z)J{xj,mj)J{xk,mk) — u'^{xk — Xj + z)J{xk,mk)J{xj ,mj) 

— u'^{xi — Xj + z)J{xi,mi)J{xj,mj) 
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Recall that we by our convention, we do not display the dependence of J and Jon the i-variable. 
To bound these terms, we require two variants of Lemma 13.31 

Lemma 3.6. There exists a collection of constants C : ^ (0, oo) such that, for any continuous 
functions t, v, ai, a2, as : M'' — > [0, oo) and any z € M*^, 



(3.51) Ejv / dt V -f{mi,mj)t 

Jo ■ .ttI ,.s 



V 



i,j,keiq{t) 

ai{xi)a2ixj)a3{xk)l{mi = ni,mfc = n^} 



"2 «J,fcG/q(0) 

where A'^^ : M^'^ [0, oo) is given by 



3d+2 

xi - z'\^ \X2 - yr , |X3 - y'" ~ 



v'l2 „,|2 „/|2\ — 2" 



^^,n2,n3(^i,^2,X3) = ^0 (3ci)7(ni , n^) y^^^ ^^^^ + + 

(3.52) t{^- ^--^) v{^- ^— ^) oi {z')a2{y)a3{y')dz' dydy, 

with 7 as in Lemma 3.3. 

Lemma 3.7. There exists a collection of constants C : ^ [0, oo) such that, for any z G M*^, any 
continuous functions v,w ^ [0, oo) and another (ai, 02, as) : [0, 00), 

(3.53) En j dt ^ -f{ni,nj)Ve{xi - Xj)v(^— — — ai(xi)a2(xfc)as(x/) 

IL{m, = ni,m, = n3,m,=n4} 

^ C'ni,n3,n4e *■ ^ ^ Ejy B^^,i^.,^^,^^{xi^Xj,Xk,Xl) 

IL{mj < ni,mj < n2,mk < n^^mi < ^4}, 
w/iere -B^^ „2,n3,n4 : I^^'' ^ [0^ o^) ^-^ ff^^en by 

(3.54) ^ni,n2,n3,n4(3;i'3;2,2:3,2:4) 

|xi — z| ^ |x2 — z I ^ |xs — 2/r \X4 — y ' ^ 



d{ni) d{n2) d{n3) ^(71,4) 



' a,' 



j{ni,n2)Ve{z - z')v(^- — ^ ~)'^{~ ^ -^ai{z)a2{y)a-i{y')dzdz'dydy' , 

where the function 7 : ^ (0, 00) satisfies 

n27(ni,n2 + nsj max|l, |^ — ^ — J j < (n2 + n3j7(ni, 712), 



The proof of Lemma 3.5 is identical to that of Lemma 3.3. The proof of Lemma 3.6 is very 
similar to the proof of Lemma 3.3 and is omitted. 
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We now bound the three terms. Of the first two, we treat only Ai, the other being bounded by 
an identical argument. By multiplying out the brackets appearing in the definition of Ai and using 
the facts that 

^^(3;, - X,- + z) = e2-'^n(^i^^i±^) , 

and 

V,,n^(x. - X, +z) = ei-'^V^x(^i^^^ 
we obtain that Ai is bounded above by 

(3.55) Ce2('^-3) ^ |Vn|(^^^^^)|Vn|(^l^^^)j2(x„m,) 

i,j,k£lq 

\j{xj,mj)\\j{xk,mk)\TL{mi = Mi,mj = ruk = M2} 
+ C62('^-2) ^ ^l^ x,-x^, + z yi^ x,-X!, + z -^ |VJ(x„ m,)P 

i,j,kelq 

\j{xj,mj)\\j{xk,mk)\TL{mi = Mi,mj = = M2}. 

Let us assume that z = because this will not affect out arguments. We are required to bound the 
quantity appearing in the statement of Lemma 13.61 for each of the following cases: 

(3.56) (t, V, ai, 02, 03) G | (|Vu|, |Vu|, J^, | J|, | j|) , [u, u, |VJp, |</|, |^ 

Recall that each of the test functions J, J, and their gradients, is assumed to be uniformly bounded 
with compact support. To each of the two cases, Lemma l3.6l applies. For either of them, the right- 
hand-side of the inequality in Lemma 13.61 may be written as a finite sum of the expectations 
appearing there, with the sum being taken over triples of given masses ni,n2 and n^. Such an 
expectation is bounded above by 

(3.57) C.3(-.J:// ^ + ^+^-"1 



„X3 V d{n2) dins) 

ii'2 

— ^ — j v\y — ^ — j /i„i (xi)/i„2 {x2)hn3 {x3)dz'dydy'dxidx2dxs, 



where K = {x : \x\ < £} QM.'^ is chosen to contain the support of J and J. We bound (|3.57j) for 
the first of the two cases in 1)3. 56|) . the other three being similar. In this case, t = v = |Vu|, the 
expression (|3.57j) equals 



Cg3(2-d) 



\Vu\ (^) \Vu\ (^-A)g(z', y, y')dz'dydy' 



where G{z',y,y') is defined to be 

v'|2 U „,|2 U ,/|2 



^/(^k^ + ^(^ + ^l(^) ' h^n.i^l)hn,i^2)KAx3)dx,dX2dxs. 

As in the proof of the fifth part of Lemma 3.5 we can show that G G Lf^^. To bound the expression 

\Vu\ (^^) I Vu| (^^) dz'dydy', 
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note that, for z' £ K, 



dydy' I Vu I 



z -y 



z'-y' 



< 



a\<2e 



\Vu\ I - Ida 



Using \Vu\{x) < C\x\'^-'^ yields a bound C£2e2{rf-i) for (Tmi) . Hence, the expression appearing m 
(jTFTl) is bounded above by Ce3(2-d)g2(d-i)_ 

Similar reasoning applies in the second case listed in ()3.56|) . The bound on the term in ()3.57|) in 
this case is a constant multiple of g3(2-(i)g2((i-2)_ Applying Lemma 13.61 to the expression in H3.55() . 
we find that 

EtV r A,{q{t))dt < Ce2(rf-3)g3{2~d)^2(d-l) ^^g2(d-2)^3{2-d)^2{d-2)^ 

These bounds correspond to the two cases appearing in H3.56() . Thus, 



E 



N 



T 



Ai{q{t))dt < Ce' 



d-2 



We must treat the third term, ^3. An application of the inequality 



ai + . . . + a„ ) < n{al + . . . + 



to the bound on ^3 provided in (|3.5U|) implies that 
(3.59) 

where Yi is given by 



n=l k£lg 



mi 



-u^ixi -Xk + z)J{xi,mi + mj)J{xk,mk), 



mi + mj 

and where {Yi : i G {2, . . . , 8}} denote the other seven expressions in (|3.5U|) that appear in a sum 
over /c G /g, while Yg denotes the last term in 1)3. 5U() that does not appear in this sum. There are 
nine cases to consider. The first eight are practically identical, and we treat only the fifth. Note 
that 

2 



(3.60) 



Ce^''-^ a{mi,mj)Vei 



Xj X ^ 



[E 

k,l&I„ 



U 



U 



j'^{xi,mi)J{xk,mk)J {xi , mi ) 



In the sum with indices involving A;, / G /g, we permit the possibility that these two may be equal, 
though they must be distinct from each of i and j (which of course must themselves be distinct by 
the overall convention). 

Note that the expression 1)3.60(1 appears in the statement of Lemma [3. 71 provided that the choice 

{v,w,ai,a2,a^ = [u,u,j'^, | j|, |j|^ 
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is made. Again we set z = because this does not affect tfie estimates. Given tfiat tfie support of 
the functions ai, 02, as : W^'^ — > [0, 00) are bounded, we must bound 

fSfin V /" /" / l^i-^l^ \x2-z'\'^ \x3-y\'^ \xi-y'\'^ \-^d+i - z' \ 



n2 



(z — ^\ ^ z — y \ 
j uy j (xi)/i„2 ix2)hng ix3)hn^ {x4)dzdz' dydy' dxidx2dx3dx4, 



for a compact set L. This expression is bounded above by 
(3.62) / v(^^)u(^)u(^^^)dzdz'dydy' 



LA \ e 

„Jk4\ d{ni) d{n2) ^(ns) ^(714) / 

Ki {Xl)hn2 {X2)hn-j, ix3)hn^ {Xi)dxidX2dX'idXi, 

which is less than 

C j^v(^)u[^)u{^)dzdz'dydy'. 

The proof of this follows the proof of the fifth part of Lemma 3.5; we use the elementary inequality 
abed < {a? + b'^ + + d'^Y ^"^^ the fact that the function 

k{x) = Y.d{nf^-^/^ I K{y)\x-y\^'^-''dy 

n '' 

is locally bounded. Since f^u(^^^^dy is bounded above by Ce'^"^, we have that 

V{^)u{^)u{^)dzdz'dydy' < Ce^^'-^^ j v{^)dzdz'. 

This is at most Ce^'^~'^. Applying Lemma 13.71 we find that the contribution to 

En [ A3{qit))dt 
Jo 

arising from the fifth term in (|3.59|) is at most (7g2d-4g4(2-d)g3d-4 _ 

We now treat the ninth term, as they are classified in (|3.59p . It takes the form 

g4d-8 a{mi,rnj)Vs(^Xi — Xj^u^ (^Xi — Xj + z^'^ J {^Xi,rni^'^ j(^Xj,mj)'^ . 

This is bounded above by 

Ce^"^'^ ct{mi,mj)Ve(^Xi - Xj^, 

because < Ce^~'^ by (3.20). The expected value of the integral on the interval of time [0, T] of 
this last expression is bounded above by 

rT 



Ce^'^-^En [ dt V] a{mi,mj)Vs{xi-Xj) < Ce" 
Jo 



.d-2 

where we used Lemma l3. II for the last inequality. This completes the proof of l|3.11() . 
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3.5. Using the estimates. We now use the estimates derived in the preceding sections to prove 
the Stosszahlansatz. The precise form of the Stosszahlansatz that we derive is now stated. 

Proposition 1. Recall from \2.S\) that 



(3.63) 



Q(0) = e'^ 2 ^ a{mi,mj)Ve{xi - Xj)J{xi,mi)J{xj,mj), 



where, as before, we simply write J{xi,mi) and J{xi,mi) for J{xi,mi,t) and J{xi,mi,t). Let 
T] : M.'^ —)■ [0, oo) denote a smooth function of compact support for which f^d ri{x)dx = 1. We have 
that 



(3.64) 



Q{0)it)dt 



dt / du;P{Mi,M2)J{uj,Mi)J{io,M2) 



d~2 



d-2 



j&Iq;mj=M2 



i€lq:mi=Mi 

+ Err(e,5), 



Xj — to 



where the constants /3 : — > [0, oo) were defined in 111.^^) . and where the function Err satisfies 



hmlimsupE/vrlE'rrfe, 5) I = 0. 



Proof. Note firstly that the inequahties ()3.9() . p.lOf) and (|3.11|) imply that, for large T, 



lim limsupEjv 



f Hu{t)dt+ [ Hi3{t)dt 
Jo 



That is, 



lim lim sup e^^'^^^^ETv 



(3.65) 



rT 

I dt a{mi,mj)J{xi,mi)J{xj,mj) 

|^y^(Xj - Xj + z) - V^{Xi - Xj] 
rT 

I a{mi,mj)J{xi,mi)J{xj,mj) 

Ve{Xi — Xj + z)u^{Xi — Xj + z) 



0. 
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This implies that 



j-T 

(3.66) e2('^-2) / dt V a{mi,mj)V 
Jo , f-f. . 



(xi - Xj)J{xi,mi)J{xj,mj) 



r-T 

• I dt a{mi,mj)V'^{xi — Xj + z)J{xi,'mi)J{xj,mj) 



\l + e'^-'^u'{xi-Xj + z) 



+ err(e,z) 



T 



£2(^-2) dt a{mi,mj)e-''Um,, 



dTT I ^3 + ^ 



J{xi,mi)J{xj,mj) + err(e,2;), 



where 
(3.67) 

and where err satisfies 
(3.68) 

Writing 
(3.69) 



1 + UjYii,m2 (''^) 



hm hmsupETv e?'?'(e, -2) = 0. 



Xi Xj -\- Z 



J{xi,mi)J{xj,mj), 



it follows from 1)3. 66() that 

(3.70) 
Note that 



f Q{0){t)dt= [ Q{z){t)dt + err{e,z). 
^0 



Xi - Xj + Z2- Zl 



J{xi - zi,mi)J{xj - Z2,rnj) 



Q{z2-zi) -e'^ ^ a{mi,mj)e '^Um,,', ^ 

^ d-2 ^ \ -2tt + Z2 — Zi\ 

- ^ Q(mi,mj)e ^Um„m,{^ j 

(3.71) |(J(xj - zi,mi) - J{xi,mi))j{xj — Z2,mj) + J{xi,mi){j{xj - Z2,mj) - J{xj,mj)) 

■ Xi - Xj + Z2 - Zl 



< e'* ^ ^ a{mi,mj)e ^Um„r, 



zi\ + \z2\)K{xi,mi)K{xj,mj), 



for small zi, Z2, where K, K : N x 'K'^ x [0, 00) [0, 00) are two smooth compactly supported test 
functions satisfying 

K,K < 2(||VJ|| + ||J||)l{x G M"' : d{x, suppJ U suppJ) < l}. 

Note that the final term in (|3.71|) is equal to (l^^il + \ z2\)Qk{^2 — zi), where Qk and Qj^ are defined 
by H3.63() and 1)3. 69() with K and K replacing J and J. By 1)3. 7U() applied with Qk{0) and Qk{z) 



37 



in place of Q{0) and Q{z), and Lemma 1X11 we find that 
(3.72) 

'Xi- Xj + Z2 - Zl 



rT 

/ a{mi,mj)e~^U„ 
Jo 



By dTTTl) and (HTT^ . 

(3.73) Q(z2 - Zl] 



K{xi,mi)K{xj,mj) 



< Z\\K\\ \\K\ 



Xi - Xj + Z2 - Zl 



J{xi - zi,mi)J{xj - Z2,mj) + err{e, zi, Z2), 
where we have that Eiv| err (e, zi, 2:2)! < Z\\K\ \ | |K|| (|zi| + 1^2 1) . By p.7U() and 1)3. 73(1 . 



(3.74) 



Q{0){t)dt 



rT 

I dt J{xi — zi,mi)J{xj — Z2,mj)a{mi,mj) 



Urr 



^'''-''^-^''^-''^^S-'^[^)6-S[^)dzidz2 + err{e,6) 



/ dt duJiduj2 U, 
Jo JR^d f-f 



T 



dt I duJiduj2UMiM2[— — 



a (Xi - CJl 



iGlq:mi=Mi 

where the function err satisfies 



^ ^'^ ^a{mi,mj)J{uJi,mi)J{uJ2,mj) 

+ err{e, 5) 
a(Mi, M2) J(u;i, Mi)J(cJ2, M2) 

d-2 Sr^ x-d^(^j - ^2 



j&Iq;mj=M2 



+ err (e, S) 



hm hm sup Eat I err (e, (5) I = 0, 

and where in the last equality, we made use of the fact that the test functions J and J take non-zero 
values only on particles of a given mass, respectively Mi and M2. Recalling ()1.8() and ()3.67|) . and 
making use of the continuity of these test functions, we obtain the statement in the proposition with 
an error term of the form Err(e, 6) = err(e, 5) + err'(e, 5), where the latter error err' = 
This completes the proof of Proposition ^ 



4. Uniform boundedness of the macroscopic densities 



Before discussing the main goals of this section, let us review what has been achieved so far 
and what remains to be done in order to complete the proof of the main Theorem 1.1. Our 
Stosszahlansatz , stated in Proposition ^ allows us to replace the microscopic coagulation term 
with an expression that involves {fn*'n^){fm.*''l^) with an error that goes to zero as e and 6 go to zero 
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in this order. Here is the microscopic density of particles of size n and rj^ is an approximation 
to the identity. After passing to the hmit e ^ 0, we will have expressions of the form 

(4.1) J{fn*v^)ifm*V^)J dx + Err{d), 

where /« now represents the macroscopic density, J is a bounded continuous test function, and 
Err{5) represents an error term that goes to as (5 ^ 0. Recall that by the Lebesgue's theorem, 
the sequence {fn * Vi^} converges to fn almost everywhere. To derive the PDE (1.6), we need to 
show that in (4.1) we can pass to the limit (5 — and replace /„ * 7]^ with /„. This can be achieved 
if we can establish the uniform integrability of the family 

(4.2) {(/n*r?*)(/^*r?*) :<5>0}. 
To prepare for this, let us consider the transformation 

(4.3) q{-) ^ 9n{dx,dt) = e'^-^ ^ 5^.^t){dx)l{mi{t) = n)dt G M, 

where A4 denotes the space of measures Unidx, dt) on the set x [0, T] such that /itn > for every 
n, with Hn (M*^ X [0, T]) < TZ. We use this transformation to define the probability measure Vn 
on as the pull-back of Pjy. We equip M with the topology of vague convergence and Ai^ with 
the product topology so that A4^ is a compact metrizable space. We now outline our strategy for 
the proof of the uniform integrability of (4.2). Let V denote a limit point of the sequence Vn- 

(1) We show that V is concentrated on the space of measure sequences /j, = (/i„ : n G N) such 

that /ii(dx, dt) = fi{x, t)dxdt with /i(.x, t) < ki almost everywhere, where ki = max^^ hi{x). 

(2) Assume that for every / < n, there exists a constant ki such that the measure V is concen- 
trated on the space of with iii{dx^ dt) = fi{x, t)dxdt for a function // satisfying fi{x, t) < ki 
almost everywhere for every I < n. We then deduce that fin{dx,dt) = fn{x,t)dxdt with a 
function /„ weakly satisfying 

(4.4) ^{x, t) = din)Afnix, t) + QUf){x, t) - Y^, 

where is defined as in (1.6) and Yn > 0. In other words, we can determine the form of 
the coagulation gain term but all we can say about the loss term is that Y >0. 

(3) Under the assumptions of the previous step, we show that the function /„ satisfies 



J J fn{x,t)^ exp{—\x\)dxdt < OO, 



10 

for every p > 2, P-almost surely. 
(4) Under these same assumptions, we show that with probability 1 with respect to V, the 
density /„ satisfies 

(4.5) ^{x, t) = d{n)AUx, t) + QUmx, t) - Q^(/) (x, t), 

where Q" is defined as in (1.6) and 

n 

Q2(/) = Vn{x,t) J2 /3{m,n)fm{x,t)+Xn. 
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with Xn > 0. That is, at this stage we can only determine the form of the coagulation loss 
term up to mass n and for the rest we can only assert that Xn > 0. 

(5) The equation (4.5) is good enough to deduce that there exists a constant A;„ such that 
fn{x,t) < kn, T'-almost surely. 

(6) Applying an induction on n we show that UnidXjdt) = fn{x,t)dxdt with /„ solving (1.6) 
weakly and that fn{x,t) < kn for every n, P-almost surely. 

The main objectives of this section are Steps 1, 3 and 5 of the above outline. The remaining 
Steps 2, 4 and 6 will be carried out in Section 5. 

Firstly, we establish the claim we made in Step 1. Throughout this Section, we let ip : [0, oo) — > 
[0, oo) denote a convex and increasing function with i/j{0) = 0. 

Definition 4.1. For n G N, write f^{n,x,t) = f^{n,x;q{t)) where 

f{n, x; q) = e'^-^ J^V^x - x,)]l{n, = n}, 

ieiq 

where rj : ^ [0, oo) is a smooth function of compact support with J rjdx = 1 and r]^{z) = 
5-'^r]{z/5). 

Note that Qn *x — f^{n,x,t)dt where gn was defined in (4.3). 

Lemma 4.1. Suppose that ip : [0,oo) [0, oo) is, in addition to its other properties, chosen so 
that V'" : [0, oo) — > [0, oo) exists and has compact support. Then, for any T,s € [0, oo) satisfying 
T>s, 

Em f ^{f\l,x,T))dx<EN f i^{f\l,x,s))dx + {T-s)^ f \Vv\^dx. 

Proof Prom the identity ENF{q{T)) = Ejv-F(g(s)) + Ejv EF{q{t))dt, we deduce that 
(4.6) En ^l^{f\l,x,T))dx = Em i^{f\l,x,s))dx + Ai + Y^Ai^^. 

JRd JRd 

where 

= d(l)Ejv [ dt [ y] A^X/^(l,x,t))l{mi = l}dx 

Js JRd ^ 

and where, for example, 

^1,1 = «(1) l)EAr / dt I 2, Ve{xi—Xj)Di^j{x,t)\\mi = mj = l>dx, 
•Is JRd ■ f-f ' ^ J 

where Dij{x,t) equals 

^ (/(I, X, t) - e''-^r]'ix - X,) - e^-Wx - x,)) - V (/(I, x, t)) . 

From the monotonicity of ^, we know that Dij < 0. This in turn implies that, for each n G N, we 
have that „ < 0. 
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In bounding ^i, we make use of the equality 

(4.7) ^.Mf') = V'"(/')|V,,/Y + V''(/')A.J^ 
Since ^" is of compact support, we find that, for any given t € [s,T], 

(4.8) [ En r{fHhx,t))\V,j\l,x,t)\\{mi = l}dx 



< 



CEn [ V|V,,/^(l,x,t)|'l{mi = l} 



where in the last inequality, we used the fact that the total number of particles in the model is 
bounded above by Ze^^'^. Note also that, for t S [s,^], 

(4.9) / ENy^i^'{f\l,x,t))A.,j\l,x,t)l{mi = l}dx 

EjvV' {f^h X, t)) AJ^l, x,t)dx = - [ ENip"{f^(l, X, t)) \V^f^{l, X, t)fdx < 

Applying the equality (|4.7|) to the term Ai in (|4.6j) and using the bounds (|4.8|) and H4.9j) . we deduce 
the statement of the Lemma, each of the terms {^i,n : n G N} being non-positive. □ 

We are now ready to establish the claim of the first step in the outline. 

Lemma 4.2. Let the measure V be a limit point of the sequence {"Pat}- Then the measure V 
is concentrated on the space of fj, such that fj,i{dx,dt) = fi{x,t)dxdt with fi{x,t) < ki almost 
everywhere, where ki = maxj, hi (x) . 

Proof Fix a continuous function 7 : [0, 00) [0, 00) of compact support with J jdt = 1. It is 
straightforward to show that the map 

(4.10) ly ^ F{v) := V (y r]\x - y)^{t)u{dy, dt)^ dx, 
is lower semicontinuous with respect to the vague topology. By Lemma 4.1, 



VNidfj) vHx-yh{t)l^i{dy,dt)^ dx<EN j v(/(l,a:,0)) dx + C 



§d+2 ' 



where '0 satisfies the conditions of Lemma 4.1. Choose A; G M to satisfy k > ki = max^, hi{x). By 
approximation, we may choose ■0 : M ^ [0, 00) to be ip{z) = (z — fc)-(_. With this choice, we can 
show that the right-hand-side of 1)4. llj) goes to in the limit N ^ 00. The left-hand-side being 
non-negative, it has zero limit also. Using the lower semicontinuity of (|4.1()|) we deduce 

V{dfi) [ V ( / v\x - y)-f{t)niidy, dt)] dx = 0. 
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Since the function 7 is arbitrary we deduce that ■p((i;u)-almost surely, the measure fii has a density 
/i and that fi{x,t) < ki, for almost ah {x,t) G M'^ x [0,oo). □ 

We now turn to Steps 3 and 5 of the outline. The second part of Lemma 4.3 below will take care 
of Step 5. A microscopic version of the first part of Lemma 4.3 is out Step 5 and will be carried 
out in Lemma 4.4. Indeed Lemma 4.3 is a variation of a theorem of Wrozosek j^j regarding the 
boundedness of the solutions of (1.6). We give a detailed proof of Lemma 4.3 even though the proof 
is a straightforward adaptation of the arguments in [Hj . 

Lemma 4.3. Let J : M'^ — > (0, 00) be a smooth function with f J dx < 00 such that for a constant 
R> ^, we have that |VJ(x)| < RJ{x) for every x. 

• Assume that {fj : R"^ x [0, 00) [0, 00), j G {1,2, . . . denotes a collection of functions 
satisfying (4-4) with Y nonnegative. Assume also that there exist constants {ki, . . . , kn-i} 
such that sup^ ^ fn(,x, t) < kj for j < n. Then for every p > 2, 



(4.11) 



sup / fn{x,tY J{x)dx < 00. 



In addition to the above assumptions, suppose that fn satisfies (4-5). Then 



(4.12) 



sup/„(x,t) =: kn < 00. 

x,t 



Proof We first establish ()4.11|) . Let us write ^pp{z) for z^. Note that 



^ l__i'p{fn{t,x))j{x)dx 



'^'p{fn{t,x)) 



d{n){Afn){t,x) 



n-l 



+ ^ a(m, n - m)fm{t, x)fn-m{t, x) - Y 



J{x)dx, 



We have that 
(4.13) 



/ i;'p{fn{t,x)){Afn){t,x)j{x)dx 

i;';{fn{t,x))\Vfn\\t,x)j{x)dx 
i^'p{fn{t,x))Vfn{t,x) ■ VJ{x)dx. 
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We find that 



< R 
R 

< — 



ijp{fn{t,x))Vfn{t,x) ■VJ{x)dx 



i'p{fn{t,x)) 

R f \Vfn{t,x) 



. %{fnit,x)y 

^';{fn{t,x)) 



R f 

■J{x)dx + — / 



i?2 



-1pp{fn{t,x))j{x)dx 



(4.14) < i?3 / Vp(/n(i,x))J(x)dx + — / |V/„(t,x)|V;'(/n(t,a;))J(x)dx. 



2i? 

In the first inequahty, we apphed the bound |VJ| < RJ{x). In the third, we made use of the 
inequahty 

/,/|2 



P 



yjp p — 1 
From ()4.13() and H4.14() . we learn that 

i;;{fn{t,x)){AU){t,x)j{x)dx < (-L-l) f i;;{U{t,x))\Vfn\\t,x)j{x)dx 

+R^ [ ijp{fn{t,x))j{x)dx 



(4.15) 

We find that 



< 



/ ^p{fn{t,x))j{x)dx. 



^ / 1pp{fn{t,x))j{x)dx 



n-1 



< 



< 



S2 0.(111,11-771) / ll;' {fn{t,x))f.m{t,x)fn-m{t,x)J{x)dx 

+ R^d{n) I i;p{fn{t,x))j{x)dx 
A [ i;'p{U{t,x))j{x)dx 
+R^din) [ iPp{fn{t,x))j{x)dx, 



where the uniform bounds sup /j</cj,iE{l,...,n — 1} provided by the hyptohesis, were apphed 
in the second inequahty, and where the constant A was defined by 



n-l 



A = a (m, n — m) kmkn- 



m=l 



Note that 
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implying that 



4 / i>p{fn{t,x))j{x)dx 
at jf^d 

< B [ i)p{Ut,x))J{x)dx +C, 



where 
and 

This impHes that 

sup / iljp[ fn{t,x)). J {x)dx < oo, 

tG[0,T] JR'* 

and concludes the first part of the Lemma. 



B = pA + R^d{n), 
C = pA I J{x)dx. 



We now turn to the second part of the Lemma. Note that 
d 

(4.16) —fn{x,t)=d{n)Afn{x,t) + ^a{m,n-rn)frn{x,t)fn-rn{x,t) - a{n,n)fn{x,t) - X, 

where X >0. We now choose ipp : [0,oo) — > [0,oo) according to iljp{z) = [{z — 1)"*"]^. A repetition 
of the proof of (|4.15j) yields 

/ i;;{fn{t,x)){Afn){t,x)j{x)dx < i;';{U{t, x))\V fnf [t, x) J{x)dx 

[ 'll^p{fnit,x))Jix)dx 



(4.17) < / i;p{fn{t,x))j{x)dx. 
It follows that 

(4.18) - / i^p{fn{x,t))j{x)dx 

at J^d 

< a[ Mfn{x,t))Jix)dx + din)R^ [ i)p[fn{x,t))J{x)dx 
-a{n,n) / fn{x,tf'ip'{fn{x,t))J{x)dx. 

Note that 

pP 

ijj'p{a) < rV'p(a) + 

implying that 

(4.19) 4 / ^Pp{fn{x,t))J{x)dx < {Ar + d{n)R^) [ i;p{fn{x,t)) J{x)dx 

at ./iRd ./rod 



-a(n,n) / fn{x,tftp'[fn{x,t))J{x)dx. 

jRd 

+A^^ f J(x)dx. 
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On the other hand, a'^ip'{a) > pip{a). Hence, the right-hand-side of (|4.19|) is bounded above by 
(Ar + d{n)R^) [ '4)p{fn{x,t)j{x)dx + A-^ [ J{x)dx - pa(n,n) [ 7pp(fn(x,t))j{x)dx. 

It suffices to estabhsh (|4.11j) for sufficiently large p. We choose p large enough so that our choice 
r = A"^ (^pa(^n,rij — d{n)R^) for r is positive. Such a choice yields 

^ l^p{fn{x,t))j{x)dx < —— ^ ^ I J{x)dx. 

jRd [pa[n,n) - d[n)R-^y JR'' 

From this, it follows that 

/ i^p{fn{x,T))j{x)dx < [ ipp{K{x))J{x)dx + T— -— ^— / J{x)dx. 

JRd jRd [a[n,n) - d{n)R-^)'^ JR'^ 

By taking the p-th root of this inequality and then taking the value of p to infinity, we deduce that 

(4.20) sup/„(x,t) < sup/i„(x) + —, r-. 

x,t X a[n,n) 

□ 

We now concentrate on the third step of the outline. Recall the random measure gn that was 
defined by (4.3). For our purposes we need to regard the microscopic density as a member of the 
smaller space A = L°°([0,r];MzM) C M where Mzi^"^) denotes the space of nonnegative 
measures 7 with 7(1^^^) < Z. More precisely, the measure gn{dx,dt) € A is absolutely continu- 
ous with respect to the Lebesgue measure for almost all x and if gn{dx,dt) = gnidx,t)dt, then 
<^„, (M'^,t) < Z for almost all t € [0, T] and every n. Since we have a uniform bound Z on the total 
measure of g{dx, t) for almost every t, onee can readily show that indeed A is a closed subset of the 
compact metric space Ai. Consider the transformation 

(4.21) q ^ {g{dx, t),X{dx, dt)) £ M = x , 
where the random measures ig,X) = {gn,Xn : n € N) are given by 

(4.22) gn{t,dx) = e'^-^ ^Mt)i'^'')'^i^ii*) = n)' 

Xnidx,dt) = e'^^'^ ^ a{mi{t),mj{t))Ve{xi{t)-Xj{t))l{mi{t) + mj{t) = n) 

i (mi{t)6^^(^t) + "^^(i)^,^) (dx) dt, 

and M denotes the space of nonnegative measures on M'^ x [0, T]. The transformation (|4.21|) defines 
a probability measure Pm for each on the space A4. This measure is an augmentation of the 
measure Vjy that is defined on the space M^. We equip A4 with the product toplogy with A and 
Af both equipped with the vague topology. Certainly A^ is a compact metric space but Af is only 
a separable complete metric space. As an immediate consequence of Lemma 3.1, we learn that the 
sequence {Pn} is indeed tight. More precisely, if Afi denotes the space of measures X € Af such 
that X(M'^ X [0,r]) < I, then sup^pAr(AAf) > 1 - C/l, for a constant C. 
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Lemma 4.4. Let J be as in Lemma 4-3. Let V be a limit point of the sequence {Vn} and write V 
for the g-marginal of V. Assume that there exist constants ki,k2, ■ ■ ■ , kn-i, such that the measure 
V is concentrated on the set of measure sequences g = {g^ : n G N) such that gj{dx,t) = fj{x,t)dxdt 
with sup^f fj{x, t) < kj for j = 1,2, . . . ,n — 1. Then gn{dx, dt) = fn{x, t)dxdt for V-almost surely 
and 



\ sup 



fn{x,tfJ{x)dx ) V{dg) < oo. 
J 

Moreover, we have that Xn{dx,dt) = ^n{x,t)dxdt, with 



n-l 



sup7„(x,t) < ^ /3(m,n-m), 

{x,t) 



m=l 



V -almost surely. 



Proof The statement of the Lemma asserts a microscopic analogue of the first part of Lemma 
4.3. Again we write '^p{z) for z^. Let ^ : M ^ [0, oo) be a smooth function of compact support with 
J £ = 1 and set f^{n,x,t) = f^{n,x,s)£{t — s)ds. A straightforward calculation yields, 



(4.23) f\n,x,t)< - f\n,x,0)iit)+ d{n)(y2 A,j\n, x, s))iit - s)ds 

/•oo 

+ / e*^"^ ^ Q;(m, n — m) ^ Ve[xi{s) — Xj{s)) 

(%^(.-..(.)) + %^(x -.,(.))) 
l{mi(s) = m, mj{s) = n — m}i{t — s)ds 

roo 

+ / Ms{x)e'(t-s)ds. 
Jo 

Here we have an inequality because we only used the gain term in applying the coagualtion operator 
Ac, and Mg is a martingale satisfying 

Ejv[m,(x)21 =Ejv / Ai{q{t),x)dt + EN [ A2{q{t),x)dt, 
L -I 7o Jo 

where Ai{q,x) and A2{q,x) are respectively set equal to 



Ai {q,x) = £^(^-2) ^d{n)\ V^ri^ | (x - Xj ) 1 (mj = n) 



and 



(4.24) A2[q,x)= ^ a(mi,m,)V;(x, -X,) 



' K{x — Xi,mi + rrij) H ^ — K (x — Xj , rui + mj) 



n 2 

—K{x — Xi, rrii) — K(^x — xj, rrij) 
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where K{x — y,m) = if{x — y)\{m = n). We can readily show, 

(4.25) ! A^{q,x)dx < Ce'^^'^-^^d-^'^-^Y d{n)l{mi = n) <Ce'^-^d-'^-\ 

[ [ A2{q{t),x)dxdt < Ce^('^-2)(5-2'^Ejv / V a{mi{t),mj{t))Ve{xi{t) - Xj{t)) 

Jo Jo , r 



l(mj(i) + mj{t) < n)dt 



where we have used Lemma 3.1 for the last inequality. 
As a result, 



dx < Ce'^-^S-^'^-'^. 



(4.26) En I fM,(a 

jRd. I 

We certainly have, 

/ -rMf\n,x,t))J{x)dx <Bi + B2 + B3 + B4, 

jRd dt 



where 



(4.27) Bi = -iit) [ iP'p{f{n,x,t))fin,x,0)J{x)dx 

JRd, 

B2 = d{n) f ^'j,{fin,x,t))(y2A,J\n,x,t))jix)dx, 



R — ^<^-2 

n?, — e 



poo r ^ 1 

/ ds i{t — s) / ip!p[f^{n,x,t)) ^^ a[m,n — m) 
Jo Jr'' 



(4.28) Ve[xi{s) — Xj{s))r]^ (x — Xi{s))l{^mi{s) = m, mj{s) = n — J(x)dx, 

and 



poo r 

(4.29) EnB4 < Ejv / dsi'{t-s) ip'(f{n,x,t))\Ms{x)\J{x)d 

Jo JRd 

< C(^-P"'+^ r ds f'it -s)\ [ En [Ms{x)] ^J{x)d 

Jo [JRd- 



1/2 



< C^(t)£'^/2-l^-l-pd^ 

where we used < Zb"'^ for the the second inequality. 
On the other hand, 

B2 = d{n) [ i;'p{fin,x,t))Aj\n,x,t)Jix)dx = B2i+B22, 
Jr^ 



where 



B2i = -d{n) [ ^j;{f^{n,x,t))\vj^{n,x,t)fj{x)dx 

JRd 
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and 



Note that 



(4.30) 



B 



22 



d{n) / i^'p{f{n,x,t))\/,f{n,x,t)-V,J{x)dx. 



I-B22I < Rd{n) 
Rd{n) 

Rd{n 



R'* i^p[f {n,x,t)) 



< 



i;'^{f{n,x,t))\ 



-J{x)dx 



+ 



1 



It follows that 
(4.31) 



1^21 < 1 



< 



1 

2R 

R^ [ ijp{f\n,x,t))J{x)d3 



^ I V, {f'{n, x, t)) I UHn, X, t)) J{x)dx 



B2i+R^ ^p{f{n,x,t))j{x)dx 



where in the second inequality, the fact that B21 < was used. 
Note that 

poo p 

B3< dsl{t-s) ijp{f{n,x,t))x\x,n;q{s))J{x)dx 



where 



n-l 



(4.32) X^{x,n;q) = e'^"^ ^ a(m,n - m) ^ F,(xi - Xj)l{ 



Xi — XijUi rrii = m,mj = n — m 



m=l i,j<^Iq 

— r/ [X — Xi) H -r] [x — xj^ 



n 



n 



In summary, the expression 



d 
di 



ipp{f^in,x,t))j{x)dx 



is bounded above by, 

(4.33) -e{t) [ i;'p{f{n,x,t))f^{n,x,0)J{x)dx 

+R^ [ ^p{f\n,x,t))j{x)dx 

+ / ds£{t-s) / 7pp{f\n,x,t))x\x,n;q{s))J{x)dx + M' 



with Af^ satisfying 
(4.34) 



E 



N 



Before sendinge — > in 1)4. 33(1 . recall the space A/j which was defined right before the Lemma and 
let us observe that the functional T : M x J\fi ^ [0, 00), defined by 




^' { I I V^x - y)i{t - e)g{dy, dO) ]ri\x- z)J{x)i{t - s)X{dz, ds)dx, 
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is continuous with respect to the weak topology. Indeed, if ^ g and X, then F'^ 

and y ^ y, where 



F 



(4.35) 



v\x-y)eit-e)g{dy,de) 



F{x) = 

F'\x) = J 7]\x-y)iit-9)g'idy,de) 

Y{x) = J rj\x-y)i{t-e)X{dy,de) 

J v\x-y)eit-e)x^{dy,de). 

Since both sequences {F^} and {V} are uniformly bounded, we deduce that 
lim / ij'{F''{x))Y''{x)J{x)dx = [ il^' {F{x))Y{x)J{x)dx. 



We now fix S and send e to 0. From (|4.33|) we deduce that the measure V is concentrated on the 
space of [gn,Xn) for which 

(4.36) ^ I i^'p{fi)Jdx < -m I ^p{f^)hiJdx + R' J ^Pp{f^)Jdx + J ^'p{f^)xiJdx, 



where, 

We now claim that 
(4.37) 



fn = 9n *x *t hi = hn v\ = X T]^ *t L 



n-1 



< /5(?Ti, n — mjkmkr, 



m=l 



"P-almost surely. By the Stosszahlansatz , as it is recorded in ()3.64|) . 



(4.38) lim lim sup / 

(5^0 N^oo Jw 



JR'^ 



K{x,t) 



n— 1 



Xn{dx,dt) - ^ P{m,n- m)gm *x r]\x,t)gn-m *x r]\x,t)\dxdt 



m=l 



dVN = 0, 



for any smooth function K of compact support. From this, the continuity of the integrand with 
respect to {g,X), and our assumptions on gj for j < n we deduce. 



/poo p I " 1 \ 

/ / K{x,t) Xn{dx,dt) — y (3{m,n — m\kmkn~m dxdt 



dV = 0, 



which implies that 



poo p 1 poo p 

/ / KdXn < l3{m,n -m)kmkn-m / / K{x, 

Jo JR'i ^0 Jm.'^ 



t)dxdt 
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for V- almost all realizations of the random measure X„. So, Xn = jndxdt, with 7^ bounded, 
implying H4.37() and the second part of the Lemma. Applying 1)4. 37p in (|4.36|1 . using V'p(^) ^ 
p^p{z) +p and Gronwall's inequality, we find that 

sup / 'ipp(^f^{n,x,t)^J{x)dx<c, 

tG[0,T] JR'i 

for a constant c that is independent of £. From this we can readily deduce, 

.up sup / *,(,„..,*(.,*,)J,.)<i.<oo. 

t6[0,T] S>0 Jr'' 

as required to derive the first part of the Lemma. □ 



5. Deriving the PDE 



The aim of this section is to carry out the second and fourth steps of the outline for deriving 
the system of PDE described in Section [l] We then apply an induction on the number n to 
complete the proof of Theorem 1.1. In each case, we choose a test function K : [0, 00) x M*^ ^ M 
and consider the functional 

Y{q,t) = e'^-^^K{t,Xi,mi), 



where K(t,x,'m) = K{x,t)Tl{m = n). Note that 



(5.1) Y{q{T),T) = y(g(0),0)+/ Ac{Y){q{t),t)dt 

Jo 

+ Ao{Y){q{t),t)dt+ I Yt{qit),t)dt + Mt, 
Jo Jo 



where 

(5.2) AcY{q,t) = e'^"^ ^ a(mi,mj)K(xi-Xj) 



m,: 



rrii + rrij 



K(t, Xi,mi + rrij^ -\ ^ K{t^ Xj,mi + rrij) 



rrii + rrij 



-K(t,Xi,mi) — K (t, Xj,mj^ 



and where 
(5.3) 



AoY{q,t) 
Yt{q,t) 



e'^ '^^d{mi)^,^K{t,Xi,mi), 

ieiq 

e'^-^^Kt{t,Xi,mi). 



The term Mt is a martingale satisfying 





= ^N [ 




Jo 
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where Ai{q,t) and A2{q,t) are respectively set equal to 



Ai {q, t) = £^('^-2) d{mi)\ Vx,i^ I' (t, x^, rm) , 



and 



(5.4) A2{q,t)= e2('^-2) a(mi,m,)K(xi-x,) 



-fC (t, Xj, rui + TTij) H ^ Xj,mi + nij) 



rrii + rrij 



mj + ruj 

—K(t,Xi,mi^ — K(t,Xj,mj^ 

We can readily show 

(5.5) Ai {q, t) < Ce^^'^-^) ^ d{mi)l{m^ = n) < Ce'^'^, 

A2{q,t) < Ce^^'^-^) a{mi,mj)V^{xi - Xj)l{mi + mj < n) < Ce'^'^, 

where we have used Lemma 3.1 for the last inequality. From these inequalities and Doob's inequality, 
we deduce that the martingale tends to zero uniformly, in the e J, limit. 

We rewrite the terms of (|5.ip in terms of the empirical measures introduced in (4.23). We have 
that 



(5.6) 

and that 
(5.7) 



^''(^) = E / K{T,x)gn{dx,T) 



[ AoYg{t)dt 
Jo 



d{n) / ^xK{t,x)gn{dx,t). 



The Stosszahlansatz enables us to rewrite the time-averaged action of the collision operator on 
y(g(t),t) in terms of these empirical measures. That is, by Proposition^ 



dtkcY{q{t),t) 



dt I dxK{t,x) 
JR'i 



(5.8) 



n-l 



P{m, n-m) [gm *x i]^) {x, t) [gn-m *x il^) {x, t) 



m=l 



2 J2 r) {gn *x V^) {x, t) {gr *x V^) {x, t) 



r=l 



To carry out the second step of the outline, note that by the assumption, the first n — \ marginals in 
the mass variable of any weak limit of the measures Vn have bounded densities {/m : [0, oo) x M'^ — > 
[0, oo) :mG{l,...,n — l}} with respect to Lesbesgue measure . By passing to the limit in low e 
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in (|5.8j) . we find that any weak hmit V is concentrated on the space oi g = {gm : m € N) such that, 



(5.9) 



gn [dx, T) K (T, x) dx 

gn{dx,O^K(^0,x^dx + d{n) / dt gn{dx,t^AK(t,x)dx 
■ Jo Jw 



T 



+ dt dxK{t,x) 



n-1 



^ /?(m, n - m) [gm *x t]^) {x, t) {gn-m *x V^) [x, t) - Xn,5 



,in=l 

+Err{6), 

where Xn^s ^ and the P-expectation of Err (^5) goes to zero as J | 0. We know by the third step 
in the outhne that gm{dx,t) = fmix,t)dx with fm being uniformly bounded for m < n. We find 
that 



gn{dx,T)K{T,x)dx = / gn{dx,0)K{0, 

(5.10) + / dtK{t,x) 

Jo 



xjdx + d{n) 



ds dt gn{dx,t^AK(t,x)dx 
Jo Jm.^ 



n-1 



^ (3{m,n- m)fra{x,t)f, 



■n—m\X, t) X^ 



where Xn > 0. From this we would like to conclude that the measure g„ = fndx and that the 
sequence {fm : m < n) solves ()4.5p in weak form. Since we do not know the absolute continuity of 
gn with respect to the Lebesgue measure yet, let us observe that if we set /„ = gn *x , then the 
sequence (/„ : n < m; fn) does satisfy an equation that is similar to (4.5). This equation is close 
enough to (4.5) so that we can apply Lemma 4.3(i) to deduce that 

sup sup J (^gn *x J dx < oo, 

for every p > 2. From this we deduce that indeed gn = fndx with /„ satisfying the same L^'-bounds. 

To carry out the fourth step, note that at the end of the third, we know that any weak limit of the 
measures Vn is concentrated on those measures on having densities with respect to Lebesgue 
measure {/m : [0, oo) x R'^ — > [0, oo) : m G N} with fm uniformly bounded for m < n and and that 
/„ G LP{R'^ X [0, T]) for each p > 2 and T > 0. Thus, the families { {fn*r]^) {fm*V^) ■ ^'^ x [0, oo) 
[0, oo), 5 > O} are uniformly integrable for m < n. When we pass to the limit in low e of (|5.8() in 
this case, we find that, 



(5.11) 



/„(r, x)K(T, x)dx 



T 



fn{0,x)K{0,x)dx + d{n) dt fn{x,t)AK{t,x)dx 



+ I dt I dxK{t,x) 



n-1 



.ni=l 



^ /3 (m, n - m) fm {x, t) fn-m {x, t) 

+ Err [6), 



- '^P{n,m)fn{x,t)fm{x,t)-Yn 
m=l 

where Yn > and the P-expectation of Err {6^ goes to zero as 5 J, 0. 
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By induction on n we deduce that the functions {fn '■ n ^ N} being known to be uniformly 
bounded ,P~almost surely. We then apply Lebesgue's theorem to deduce that, for each n G N, 
fn*V^ ^ f almost surely, as (5 J, 0. We find that the weak limit V of the sequence Vn is supported 
on those functions / satisfying the system in weak form. 

6. Potential theory 

The purpose of this section is twofold. Firstly, we show the existence of the function u that was 
used in the proof of the Stosszahlansatz . Secondly, we establish a bound on /? and show that if 
a ^ oo then /3 converges to macroscopic coagulation rate in the case of a hard sphere interaction. 
The limit can be expressed in terms of the Newtonian capacity of the support of the function V. 
We start with the statements of the main results of this section. Let y : M'^ — > M be a continuous 
function of compact support with V > and / V{x)dx = 1. We also write Kq for the topological 

closure of Uq where Uq = {x : V{x) ^ 0}. 

Theorem 6.1. For every /3 > 0, there exists a unique function u E C^(M'^) such that lim u{x) = 

|a;|^oo 

and 

(6.1) An - (3Vu = pV. 
Moreover — 1 < u{x) < for every x S M"^. 

We have applied this theorem in Section 3 with 3 = — ^— - — -. Also recall that d(n,m) is 

d(n) +d{m) ^ ' 

defined by 

(6.2) f3{n,m) = a{n,m) / V{x){l + u{x))dx. 

More generally, we define F : (0, oo) — s- (0, oo) by F{(3) = P V{1 + u^)dx with u = as in 
Theorem 16.11 Recall that the (Newtonian) capacity of a compact set ii' C R'* is given by 

(6.3) Cap(K) = inf | - / iWpdx : V G C^(M'^), V > 1 on a neighborhood of , 
or equivalently, 

(6.4) Cap(i^) = sup < ^i{K) : supp^ C K, co \x — y\'^~'^fj,{dy) < 1 for every x > , 

where cq = co{d) = (d — 2)^^uj^^ with cj^^ denoting the surface area of the unit sphere. We are now 
ready to state our second main result. 

Theorem 6.2. We always have F{j3) < Cap{KQ). Moreover lim F{(3) = Cap{KQ). 

13— *oo 

Proof of Theorem 16.11 

Step 1. Let J be a bounded continuous function with J > 0, / J{x)dx = oo and 
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/ J{x)\x\^-^'^dx < oo. Define 

J\x\>l 

(6.5) H = <u : n is measurable and / u^{x)J{x)dx < oo > . 
We then define !F -.H ^ Hhy 

(6.6) T{u){x)=co[ \x-y\^-^Viy)u{y)dy. 
Let us verify that J^{u) G H for n E "H. To see this, write 

(6.7) r{x) = co[ \x-y\^-''V{y)dy. 



As in (|3.29j) we can readily show 
(6.8) < r(3;) < cimin(|xp-'^,l) 

for a constant ci. Also, we may use Holder's inequality to assert 



n 2 



r{x) / colx-yl'-^'ViyHy) 



r(x) 

2-^^T//„.^„.2/ 



(6.9) (-^(n)(x)) 

< r(x) / colx-2/|"-'*y(y)u"(2/)a!y. 
From this and (|6.8() we deduce 

/ {T{u){x)fj{x)dx <CQCi I V{y)u^{y) I min(|xp-'^, l)[x - yp-'^J(x)(ix 
If V{y) / then jy| < Rq for a suitable Rq. If 



I{y) = I min(|xp '^,l)\x-yf '^J{x)dx, 



then we have sup I{y) < oo because 

\y\<Ro 



I{y) < + [ min([xp-'^,l)|x -yp-'^J(a;)dx 

J\x\<2Ro J\x\>2Ro 



< C \x - yl'-'^dx + C \x\'-''\x\^-''J{x)dx < oo 

J\x-y\<2Ro J\x\>2Ro 

by our assumption on J. As a result, 

{T{u){x)fj{x)dx <C [ V{y)u^{y)dy < oo 

because V is of compact support. This shows that T{u) E 7i whenever u G H. 
Step 2. Observe that 7^ is a Hilbert space with respect to the inner product 



{u,v) = / u{x)v{x)J{x)dx. 

Note that if u solves (|6.1|) , then 

n + /3(-A)'i((l + n)y) = 0. 
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As a result, 

(6.10) {id + (3T){u) = g 

where g{x) = — /3r(x) with T as in ()6.7() and id means the identity transformation. The formulation 
(|6.1U|) allows us to use some standard arguments to establish the existence of a solution of (|6.1jl . 
Note that our assumption on J implies that T G H because of ()6.8p . Hence as a strategy for 
the existence, we show that the range of the operator id + PJ^ is the space Ti. By the Fredholm 
Alternative Theorem, it suffices to show that the operator is compact and that the operator 
id + has a trivial kernel. 

Step 3. In this step we verify the injectivity of the operator id + (3J-. Let us assume that for 
some u (z ?{, (3J-U = —u. From this we deduce Au = (3uV weakly by standard arguments. Using 
Sobolev's inequalities and a bootstrap we can readily show that n € because V is continuous. 
As a result. Am = f3uV in the classical sense. We also have 

(6.11) u{x) = -co/3 / u{y)\x - y\^-''V{y)dy. 
We now assert 

(6.12) lim u{x)\xf-'^ = -cq(3 / V{y)u{y)dy, 



(6.13) lim |Vn(x)||xl^"^ < co((i-2)/3 / V {y)\u{y)\dy . 

The proof of (|6.12|) follows from 

lim \xf-'^\x - y\^-'^ = 1 

\x\~-*oo 

for every y, and the elementary inequality 

for |x| > 2i?o (because \y\ < Rq and \x — y\ > \x\/2 whenever |x| > 2Rq). The proof of (|6.13|) is 
similar. We now choose R > Rq and use Au = (3Vu to write 



uAudx = P Vu dx. 

\x\<R J\x\<R 



After an integration by parts we obtain 



|Vn| dx + vSJu ■ nda = P Vu dx, 

\x\<R J\x\=R J\x\<R 

X _____ 

where n = - — - is the normal vector and da is the surface measure on |x| = R. From 1)6. 12() - (|6. 13(1 

\x\ 

we can readily deduce 

/ uVu-nda = 0{R^-'^). 

J\x\=R 

As a result, 

- / iVupdx = P [ Vu^dx. 
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From this we deduce / \\7u\'^dx = / Vu^dx = 0. This in turn imphes u = 0. 

Step 4. In this step we verify the compactness of the operator J^. Set (/)(x) = co|x[^~'^ so that 
= (Vu) * 0. We then have VJ^(tt) = (Vu) * V(j). Let Br denote the ball {x : \x\ < R}. We 
then use the Young's inequality to write 

(6.14) \\VT{u)h2(^Bn) < \\y^^\\LHBn)\\'^4>\\mBn) < CR\\Vuh2(^Bn), 



(6.15) mu)\\LHBn) < \\Vu\\mBjmmB^) < CR^Vuh2(^B^^, 

for some constant C. Now if {un} is bounded in Ti, then by positivity of J we deduce that 

sup WunWi^^Bn) < This and (|6.14|) - H6.15() imply 

n 

(6-16) sup(||J^(u„)||i2(B^) + ||VJ^(n„)|Ii2(B^)) < oo. 



From this and Rellich's theorem we learn that {^(u„)} has a convergent subsequence in 
As in Step 1 we have 

{T{Un){x)fj{x)dx < C 



loc 



'\x\>l 

for / > 2i?o- As a result, 



x\^~^'^J{x)dx 



x\>l 



yiy)ul{y)dy, 



lim sup 



{T{un){x)fj{x)dx = 0. 



\x\>l 



This and the precompactness of {^(n„)} in Lf^J^'^) imply the precompactness of {T{un)} in Ti. 
Step 5. So far we have shown the existence of a unique solution u ^ 7i oi u + f3J-{u) = g. As 
in Step 3 we can readily show that u is a weak solution of (|6.1|) . Again we can use the Sobolev's 
inequality and a bootstrap to show that in fact u G C^(M'^) and satisfies (|f).l|) in the classical sense. 
Step 6. It remains to show that for the unique solution u of ()6.1() . we have — 1 < u{x) < for 
every x. First take a smooth function 995 : M — > (— oo, 0] such that v?^ > and 



r > -1, 

1 + r r < -1 



(5. 



We then have 
(6.17) 



Lp'^{u)\V uf' dx 



ips{u)Audx = (3 / V{1 + u)ips{u)dx 



by an integration by parts. The formula 1)6. 17() is established with the aid of ()6.12|) as in Step 3. 
Since the left-hand side of (|6.17|) is negative and {1 + u)ips{u) > we deduce 

(f'^{u)\Vufdx = / V{1 + u)(fs{u)dx = 0. 



We now send 5 — > to deduce 



0=/ \Vu\^l{l + u<0)dx= V{l + ufl{l + u<0)dx. 

As a result, on the set A = {x : 1 + u{x) < 0} we have Vu = 0. Hence u is constant on each 
component B of A. But this constant can only be —1 because on the boundary of A we have 
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1 + n = 0. This is impossible unless A is empty and we deduce that u > —1 everywhere. Finally, 
from 

u{x) = CO / V{y){l + u{y))\x - y\^'^dy, 
and 1 + u > we deduce that u < 0. □ 

Proof of Theorem 16. 2I Let us write for the unique solution of (|6.1|) . Note that by Theorem l6.H 

(6.18) -1 < li^(x) < 
for every x. Since Au^ = /?(1 + u^)V, we also have 

(6.19) u^(x) = -Co / |x -?/p-V^(dy). 
where ^^{dy) = (5{1 + u^{y))V{y)dy. From this, (|6.18|) and (|6.4|) we learn 

(6.20) Fi(3) = [ (3(1 + ul'{y))V{y)dy = /(M"^) < Cap(i^o). 

It remains to show lim F{[3) = Cap(i^o)- By (|6.2U|) the family {n^} is precompact with respect 
to the weak topology. Let be a limit point and put 

(6.21) u{x) = -co f \x- y\^-'^fi{dy). 

jRd 

Assume lim fj,^ = fi where {/U^'} is a subsequence of {/U^}. It is not hard to see 

/3'—>oo 

lim sup (x) < u{x) 

by replacing the integrand in (|6.21j) and (|6.19j) with a bounded cutoff function and passing to the 
limit. Since > —1, we deduce that u> —1. Also 1)6.21(1 implies that -u < 0. Hence — 1 < u < 0. 

In fact it is well known that lim sup (x) = u{x) except on a set of zero capacity. (See for 

example page 161 or the proof of Theorem 5.7.1 in [T].) Since Au^^ = Au = and ji^ ji, we 
can readily deduce that ^ u weakly. As a result, 

lim / {l + u'^)Vdx= / {l + u)Vdx. 
Since (3 (1 + u^)Vdx is bounded in (3, we deduce 

/ (1 + u)Vdx = 0. 

jRd 

As a result, u = —1 almost everywhere in Uq = {x : V{x) ^ 0}. Since u{x) = limsupn(y) except 

for a set of zero capacity, we learn u = —1 in Uq except for a set of zero capacity. (See, for example. 
Theorem 3.7.5 of 1 .) Let v denote the equilibrium measure of Kq. That is, supp(i^) C Kq, 

i'{Kq) = Cap{KQ), and that the function v{x) = cq \x — y\^~'^v{dy) satisfies v = 1 in Kq except 
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for a set of zero capacity. Also, since u is finite everywhere, we learn that ^{A) = whenever 
Cap{A) = 0. The same is also true for v. (See for example Theorem 5.1.9 of [1].) As a result, 

(6.22) KKo)= [ vdfi = co[ \x-y\^-'^i^idy)fiidx) = - [ u{y)v{dy). 

But ^{A) = if Cap(j4) = and —u = 1 on Kq except for a set of zero capacity. This and 1)6. 22() 
imply 

/x(Ko) = i^{Ko) = Cap(i^o). 

□ 
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